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A  Differential  Method  of  Optical  Calculation 
Part  I :  Variables 
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Elmo  Optical  Company,  Mixuho-ku,  Nagoya 
(Received  March  30,  1961) 


Abstract 

In  an  optical  system  which  has  rotational  symmetry  about  the  optical  axis, 
it  is  generally  possible  to  define  the  incident  ray  with  three  proper  quantities 
and  to  classify  various  aberration  terms  with  them. 

By  adopting  fundamental  variables  as  such  quantities  the  image  formation  of 
an  optical  system  is  discussed  theoretically  and  formulas  for  calculation  are 
derived. 

Then  subsidiary  variables  that  are  suited  to  express  the  geometrical  rela¬ 
tion  between  a  ray  and  the  optical  axis  are  defined  to  be  used  for  later  calcula¬ 
tions  as  auxiliary  quantities  especially  suitable  for  differential  calculations. 

Further  by  introducing  principal  variables  which  are  used  to  fix  the  position 
of  a  point  on  the  ray.  the  object  and  image  points  adopted  in  this  treatise  are 
defined,  and  the  quantities  corresponding  to  the  magnification  and  the  focal 
length  for  the  paraxial  ray  are  given. 

It  is  shown  that  principal  variables  in  every  space  are  connected  by  linear 
relations  and  therefore  can  be  expressed  in  differential  forms  of  the  1st 
order.  To  make  use  of  these  relations,  a  standard  optical  system  in  which  rays 
are  fixed  by  taking  an  object  point  at  infinity  and  another  point  at  the  entrance 
pupil  is  introduced. 

A  general  optical  system  in  which  an  object  point  is  located  at  a  finite  dis¬ 
tance  can  be  transformed  to  a  standard  optical  system  by  placing  an  ideal 
optical  system  in  front  of  it,  and  as  a  result  the  treatment  becomes  unified. 

In  the  standard  optical  system,  subsidiary  variables  are  more  convenient 
to  specify  the  incident  ray  and  to  classify  aberration  terms  than  funda¬ 
mental  variables. 

Besides,  formulas  expressing  relations  among  the  above  mentioned  variables 
are  derived  and  used  extensively  in  later  calculation. 


1.  Introdaction 

Basic  properties  desirable  for  variables  used  in  the  following  calculations  are 
considered.  A  Cartesian  co-ordinate  system  x,  y,  z  of  which  the  origin  is  any  point 
on  the  optical  axis  and  the  d;-axis  coincide  with  the  optical  axis  is  ad(^ed. 
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With  these  co-ordinates  two  points  on  the  ray  are  expressed  as  Q(x,y,z)  and 
Q*(x*,  Vt,  2*). 

Hereafter,  these  coordinates  and  all  the  other  quantities  concerning  the  length 
are,  if  not  mentioned,  expressed  in  terms  of  optical  path  lengths,  that  is  the  length  - 
multiplied  by  the  refractive  index. 

Let  us  call  the  following  three  quantities  that  are  given  by  the  co-ordinates  of 
Q  and  Q*  as  the  fundamental  variables 

a=y^+z\  b=yy0+zzt,  c=y**-f2**  (I-l) 

In  an  optical  system  that  has  a  rotational  symmetry  about  the  optical  axis,  the 
aberration  terms  are  classified  according  to  these  three  independent  variables  in 
the  object  space*\  In  other  words,  the  quantity  X'  that  is  suited  to  describe  the 
condition  of  the  image  can  be  expanded  in  terms  of  a,  b  and  c  with  the  use  of  the 
differential  coefficients  of  )C  by  them,  as 

+  *•*)+•••  •  (I“2) 

If  we  differentiate  X',  for  example,  by  a  and  put  fl=6=c=0,  then  </X7</a=X'a. 
That  is  to  say  the  aberration  coefficient  which  is  equivalent  to  the  expantion  coeffi- 
'cient  of  X'  can  be  obtained  with  such  differentiation. 

If  variables  such  asy,z  and  so  on  are  taken,  the  coefficient  such  as  X'a  can  not 
be  obtained  with  differentiation  of  the  1st  order  but  with  that  of  the  2nd  order.*’*’ 
Therefore  such  quantities  as  fundamental  variables  defined  in  (I-l)  are  con 
venient  to  discuss  the  aberration  of  an  optical  system  and  are  being  used  mainly 
in  the  aberration  theory'. 

To  perform  the  ray  tracing  and  the  differential  calculation  in  an  unihied  way, 
three  subsidiary  variables  closely .  related  to  the  fundamental  variables  and  more 
convenient  for  these  calculations  are  introduced  in  this  treatise. 

But  a  point  on  a  ray  can  not  be  determined  completely  with  three  variables 
alone.  For  example,  x  and  x*  are  inditerminate  even  when  the  fundamental  vari¬ 
ables  are  given. 

Hence  in  this  treatise  principal  variables  that  are  suited  to  represent  a  point 
on  the  ray  are  employed. 

Subsidiary  and  principal  variables  will  be  explained  in  the  following  two 
sections. 

1)  Masazo  Kiuti:  Kogaku  (Iwanami  Butsurigaku  Koza,  Iwanami,  Tokyo,  1941). 

2)  A.  E.  Conrady :  Applied  Optics  and  Optical  Design  (Oxford  Univ.  press,  London,  1929). 

3)  D.P.  Feder:  J.  Opt.  Soc.  Amer.  41  (1951)  630. 
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2.  Subsidiary  Variables 

To  express  the  position  and  direction  of  rays,  a  rectangular  co-ordinate  system 
^,7,  C  is  used  also. 


2*- 


Fig.  I-l 


The  origin  P  of  f ,  ij,  C  co-ordinates  is  the  point  of  intersection  between  a  ray 
under  consideration  which  is  taken  as  ^-axis  arid  the  perpendicular  from  O  to  $■ 
axis  which  is  taken  as  7-axis. 

Let  h  be  the  length  of  OP  and  let  the  direction  cosines  of  the  axes  7,  C 
with  respect  to  the  co-ordinate  system  x,y,z  be  {gn.  On,  Qis),  On,  Oa),  (OsuOMtOn) 
respectively,  and  define  A,B,C,D  and  H  by  the  following  expressions  with  the 
use  of  these  quantities 

A=gn,  B=—hgn.  C=A(/ji,  H=h*,  D=l—A*.  (1-3) 

As  is  seen  from  these  expressions,  both  A  and  D  represent  the  relation  between 
the  direction  of  the  ray  and  x-axis.  —B  which  is  hgn,  is  the  x  compcment  of  OP 
in  Fig.  I-l.  C  is  a  variable  peculiar  to  the  skew  ray  and  is  equal  to  zero  when 
gi\—0. 

It  will  be  shown  later  that  C  becomes  invariant  throughout  the  optical  system. 
//  is  a  variable  which  depends  only  on  the  length. 

We  call  these  the  subsidiary  variables,  by  which  the  condition  of  a  ray  in 
relation  to  the  optical  axis  is  determined.  There  remains  however,  a  freedom  of 
rotation  of  the  ray  about  the  optical  axis.  Hence  the  condition  with  respect  to  y- 
and  z-  directions  will  be  treated  in  the  next  section. 

H,  B  and  D  are  small  quantities  of  the  2nd  order  and  are  generally  taken  as 
three  independent  variables  corresponding  to  the  fundamental  variables  a,b  and  c. 

The  following  relation  holds  from  (1-3)  and  will  be  used  for  check  of 
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calculation. 

HD=B^+0.  (1-4) 

3.  Principal  Variables 

To  express  a  point  on  a  ray,  the  following  principal  variables  are  used  with 
which  the  object  point  and  the  image  point  of  the  ray  in  this  theory  are  defined. 


Fig.  1-2 

(insider  a  ray  before  and  after  passing  a  spherical  surface  in  Fig.  1-2.  Let 
Q{x,y,z)  be  a  point  on  the  incident  ray  and  Q'(.x',y',z')  be  the  point  of  the  inter¬ 
section  between  the  refracted  ray  and  OQ. 

If  we  put  EIE'=bOQ'IOQ,  we  obtain  EIE'=y'ly=z'lz.  This  dedgnation  of 
Q,  Q'  and  E,  E'  is  extended  to  the  whole  optical  system. 

The  plane  which  contains  jr-axis  and  OQ  (and  therefore  OQ')  is  invariant 
throughout  the  optical  system.  We  call  this  plane  the  fundamental  plane.  Q  and  Q' 
are  the  intersections  between  the  ray  and  this  fundamental  plane.  By  these  points 
which  are  the  generallized  so-called  diapoints  of  Herzbergar^\  the  object  and  image 
points  of  this  surface  are  defined.  In  the  case  in  which  the  ray  is  always  in  one 
plane,  such  as  meridional  plane,  thfey  correspond  to  the  image  points  of  the  sagit¬ 
tal  rays. 

Using  these  relations,  we  express  any  point  on  the  ray  with  the  principal 
variables  E,G  or  F  and  the  principal  variable  in  object  space  Eo,  as 

x=—GIE,  y=EolE’yo,  z=EolE-Zt,  s=—FIE.  (1-5) 

The  suffix  o  signifies  the  quantity  in  object  space  hereafter. 

In  the  object  space,  (1-5)  becomes  by  taking  Eo  and  Go  instead  of  E  and  G 
respectively  as 

x=-G^Eo,  v=yo.  s=-FolEo.  (I-€) 

4)  M.  Helzberger:  Modern  Geometrical  Optics  (Interscience  Publishers,  Inc.,  New  York, 

1958). 
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If  the  intersection  between  the  object  plane  and  the  optical  axis  is  taken  as 
the  origin,  x  and  therefore  G,  become  zero  and  only  £o  remains. 

When  the  object  plane  is  at  infinity,  becomes  to  zero  and  only  Go  remains 
as  seen  from  (1-6).  In  this  case  the  following  expressions  containing  Go  instead  of 
Et  are  used. 

x=—GIE,  y=GtlE-yo,  z=GoJE-zo,  s=—FIE  (1-7) 

where  Vo  and  Zo  become  quantities  expressing  the  direction  of  the  ray. 


4.  Linear  Relations  among  Principal  Variables 


On  the  basis  of  the  expressions  mentioned  in  the  previous  section,  it  will  be 
shown  in  this  section  that  linear  relations  hold  between  the  principal  variables  in 
one  space  and  those  in  the  other  space. 

Here  a  point  on  the  ray  for  which  the  subsidiary  variables  are  given  is  con¬ 
sidered. 

The  principal  variables  by  themselves  have  no  physical  meanings  and  only  the 
ratios  of  one  to  another  can  express  practical  quantities  as  seen  from  (1-5)  and 
(1-7).  In  object  space  Go/Eo  in  (1-6)  corresponds  to  this  quantity. 

The  expressions  of  (1-5)  are  again  considered.  We  may  express  x  as 


G  /  E 

x=  — where  E^Eo  and  x  are  unaltered  when  the  ratio  GdEo  is  a  constant, 

because  in  this  case,  from  (1-6),  x,y,z  and  therefore  the  point  renuuns  unchanged. 
Accordingly,  GlEo  similarly  to  E/Eo  remains  unaltered.  This  relation  holds  only 
when  E  and  G  are  both  expressed  by  linear  combinations  of  Eq  and  Go. 

These  linearlities  of  principal  variables  make  our  calculation  rational  and 
especially  the  differential  calculation  simple.  It  is  mainly  because  of  this  linearity 
that  the  principal  variables  are  employed  in  our  calculation. 


5.  A  Standard  Optical  System 

In  this  section,  by  assuming  the  linear  relations  which  hold  among  principal 
variables  as  shown  in  the  previous  section,  an  optical  system  named  the  standard 
optical  system  is  derived. 

It  means,  for  example,  that  the  values  of  any  two  points  on  a  ray  determines 
the  value  of  any  other  point  on  the  same  ray  by  these  linear  relations.  Take  the 
center  of  the  entrance  pupil  as  the  origin  and  consider  the  two  points 

Go=l,  Eq=0  , 

and  • 

G,=0,  £0=1  (1-8) 
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on  the  incident  ray.  Of  th^,  one  is  the  object  point  at  infinity  and  the  other  is 
the  point  on  the  entrance  pupil.  We  call  such  an  optical  system  involving  these 
two  points  the  standard  optical  system. 

If  X  is  the  principal  variable  in  any  space,  the  value  of  X  can  then  be  represented 
by  Xa  and  Xe,  the  derivatives  of  X  by  Eg  and  Go  because  of  the  linearity  (Here¬ 
after,  the  suffices,  if  not  otherwise  mentioned,  signify  quantities  in  the  object  space). 
As  regards  any  point  in  general  with  principal  variables  Eg  and  Go  in  object  space, 
the  value  of  X  can  easily  be  obtained  as 

X=XEEg-\-XaGg  .  (I“9) 

It  will  be  clear  that  the  differential  relations  among  the  principal  variables  are 
none  other  than  what  is  described  above,  and  that  differential  coefficients  of  the 
order  higher  the  2nd  do  not  exist. 

In  object  space,  taking  (1-9)  for  X=^G  and  X=E  and  using  the  values  of  Eg 
and  Go  in  (1-8),  we  obtain  the  evidient  expressions 

Gc=l,  Ge—0 
and 

.  .  Ea=0,  Ee=1.  (I-IO) 

Let  us  employ  such  a  differential  expression  of  principal  variables  in  every 
case  hereafter. 

Thus  (1-5)  and  (1-7)  are  replaced  by  the  following  expressions 

EgX=—Go,  EaV^Vg,  EaZ=Zg,  EaS=—Fa 
EeX^=—Ge,  ^ey*  =  l/*o,  £^2*=2*o,  EeS^=—Fe  (I-ll) 

where  suflix  *  represents  a  quantity  referred  to  the  pupil. 

Then  from  (I-l)  and  (I-ll)  the  expression  of  fundamental  variables  is 

E*aO=ag,  EaEEb=bg,  E*EC=Cg.  (1-12) 

6.  Incident  Rays  and  Subsidiary  Variables 

As  stated  above,  the  object  plane  and  the  entrance  pupil  can  be  determined  by 
principal  variables.  We  now  consider  how  the  position  of  an  object  point  and  the 
direction  of  a  ray  are  expressed  by  subsidiary  variables. 

Take  a  standard  optical  system  in  which  the  object  point  at  infinity  is  on  the 
meridional  plane  and  the  center  of  the  entrance  pupil  is  the  origin  O,  and  investigate 
the  relations  between  the  rays  from  the  object  point  and  subsidiary  variables. 
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In  Fig.  1-3,  Pup2,h  and  pt  are  points  equally  distant  from  O  and  are  on  the 
plane  which  includes  O  and  is  perpendicular  to  the  principal  ray.  From  the  ex¬ 
pression  (1-3)  it  will  be  clear  that 

1  Since  the  rays  from  the  object  point  are  parallel  to  one  another,  pu  and 
therefore  Do  are  both  constants. 

2  As  regards  the  meridional  rays  onto  P2P4.  ffsi  and  therefore  C  are  both 
equal  to  zero. 

3  As  regards  the  sagittal  rays  onto  PiJPs,  921  and  therefore  Bq  are  both 

equal  to  zero.  _ 

4  As  regards  the  zonal  ray  onto  PiPzPoP*,  H  is  constant. 

As  shown  above,  the  subsidiary  variables  in  object  space,  namely  Ho,  Bo,  Do  or 
C,  are  the  variables  that  suit  to  specify  the  incident  ray  of  the  standard  optical 
system  and  consequently  to  classify  the  aberration  terms  in  this  case. 

But  they  are  insufficient  to  describe  the  incident  ray  simply  in  a  general  optical 
system.  Use  of  other  variables  such  as  the  fundamental  variables  a,  b  and  c  is  also 
effective  as  is  seen  from  (I-l).  Nevertheless,  we  shall  find  the  way  how  a  general 
optical  system  can  be  transformed  into  a  standard  optical  system  in  Section  1-9. 

7.  Fundamoital  Relations 

Fundamental  relations  between  the  quantities  referred  to  two  points  on  a  ray 
are  investigated.  •  •  •  ' 


8 


Sadao  KatO 


Z 


Fig.  1-4 


In  Fig.  1-4,  let  the  two  points  be  Q{x,y,z)  and  i/*,  z*).  Put  PQ=s, 

PQt=s^  and  Q^Q=t.  Then  there  exist  the  following  relations  between  these 
quantities  and  fundamental  variables  a,b  and  c. 

/s = (or* + o) — (XX*  +  6) 
fs* =(xx* +6)— (x**-f-c) 

/*=(x*+a)-2(xx*-l-6)-|-(x**+c)  (1-13) 

With  the  use  of  these  quantities,  the  quantities  regarding  h  and  the  direction 
cosines  are  expressed  as 

fflfii=x— X*,  fA(;2i=sx*— s*x 
tgn=v—y*,  thgi2=sy^-s^y 
tgu=z-z*,  thg2s=sz^—s*z 
thgii=yzt—y^^z 
thgsz=zx^—z^x 

thgn=xy^-x^y  (1-14) 

/*/i*=(x*+fl)  (x**-t-c)— (xx*+6)* . 

From  (1-14)  and  (1-3),  we  have  t{gi2Z—gi3.y)=yz^—y^z=thg3i=Ct,  namely 

C=gi2Z-guy.  (1-15) 

Combining  this  relation  with  (1-14)  and  using  (I-ll)  we  obtain 

CGq = /((gsiVo + 0is2o) 

CFa = hig22yo—gt»Zo) .  (1-16) 

Further  a  fundamental  relation,  in  other  words  a  quantity  invariant  throughout 
the  optical  system,  is  derived. 


A  Differential  Method  of  Optical  Calculation.  I 


9 


Fig.  1-5  shows  the  relations  of  Q  and  to  O'  and  O'*  that  are  the  image 
points  of  Q  and  0*  respectively.  The  area  of  ^OQQ^  is  t-h=OQ*.  OQ  sin  d  and 
that  of  hOQ'Q'*  is  t'-h'=6Q'.  OQ'*  sin  6. 

Taking  the  ratio  between  the  two  and  assuming  the  relation  h=h'  that  will  be 
explained  later,  we  have  E'aE'  Et'=EaEKt. 

If  this  treatment  is  extended  to  the  object  space 

EaEEt=to  (H7) 

is  obtained,  where  to  represents  a  quantity  invariant  and  its  value  will  be  deter¬ 
mined  in  the  following -section. 


8.  Relation  among  Variables 

On  the  basis  of  the  relations  derived  in  the  previous  section,  various  funda¬ 
mental  relations  among  principal  and  subsidiary  variables  will  be  obtained. 

First,  multiplying  the  first  formula  of  (1-14)  by  EaEg,  we  have  —toA=GaEB— 
GeEq.  From  (I-IO)  it  becomes  in  object  space,  as 

/o.4o=-l.  (1-18) 

Therefore,  the  fundamental  relation 

GqEe — GEEa=AIAo  (1-19) 

is  obtained.  Similarly,  by  multiplying  (1-13)  by  Ee'Eo*  and  using  the  expressions 
of  (I-l)  and  (1-14),  we  obtain 

t*oD=  E^eOo — 2£*£o6o + E^aCo 
— fioD=GEEEao — {GaEE~\'GEEG)bo-\-GaEaCo 
/*#(// — C*)=G*£Oo — 2GoG*6o4"0*oCo 
PoC*  =  OtCo—t^f 


(1-20) 
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If  (I-IO)  is  again  used  the  above  formulas  in  object  space  are  reduced  to 

fo2Z)o=Oo.  VBo=6o.  to\Ho-0)=Co,  toAo^-1.  (1-21) 

From  (1-20)  and  (1-21),  the  following  fundamental  relation  formulas  among 
principal  and  subsidiary  variables  are  obtained 

D=E*eDo-^2EeEgBo-\-E^(KHo — C*) , 

B=  EeGEDo+iEEGo-i- EoGe)Bo-\- EaGa{Ho — C*)  , 
H-C*=G'EDo+2GEGaBo+G*G{Ho-C*).  (1-22) 

It  is  seen  that  subsidiary  variables  can  be  computed  from  principal  variables 
in  any  space.  This  corresponds  to  the  fact  that,  when  the- positions  of  any  two  points 
on  a  ray  are  given,  the  relation  between  the  ray  and  the  optical  axis  is  determined. 

Therefore,  the  complete  expression  of  the  emerging  ray  in  image  space  is 
possible  with  four  principal  variables  Go,  Ea,  Ge  and  Ee.  Hence  we  can  take  these 
principal  variables  as  the  main  quantities  in  calculation  hereafter. 

Also  (1-22)  makes  the  differential  calculation  of  subsidiary  variables  simple, 
especially  that  of  the  higher  order,  because  of  the  linear  relations  among  subsidiary 
variables  D,B  and  H. 

Besides,  relations  which  are  to  be  useful  later  on,  are  here  derived.  The  com¬ 
bination  of  (1-16)  with  (1-3)  gives 

Gg=FoA-\- EfiB,  F*a~^E^GH=G*a~\~D/A\ .  (1-23) 

These  relations  are  convenient  for  checking  of  calculation.  Modifiing  (1-23),  we 
have  EGA=GG{A*-i- D)—EgB,  then  Dg  is  defined  as 

Pa=Fa — GgA=~^(GgD — EgB).  (1-24) 

For  an  axial  ray,  Bo=Do=0  and  therefore  P=0,  hence 

Fa = GgA  (1-25) 

is  given. 

Further,  from  (1-22)  and  (  -24),  GoZ>— £oB=(Go£i:— £^oGs)  (£Gfio+^e^o)= 

--^{EgBo+ EeDq)  is  derived,  then  Pa  becomes 
Ao 

AoPq^EgBq-^-EeDo  •  (1-26) 

If  Eo  is  used  instead  of  Go,  Pe  is  given  by 

-AoPe=EeBo+Eg(Ho-0).  (1-27) 

These  relations  (1-26)  and  (1-27)  are  used  later  for  ray  tracing  calculation  and 
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for  calculation  of  wave  front  aberration. 

9.  An  Ideal  Optical  System 

To  simplify  our  calculation,  a  proper  optical  system,  named  the  standard  optical 
system,  with  which  a  general  optical  system  is  unified  to  the  standard  optical 
system,  is  introduced. 


Fig.  1-6 


In  Fig.  1-6,  let  the  intersection  between  the  object  plane  OQ  and  the  optical 
axis  be  the  origin  O.  Let  Q  be  the  object  point  and  Q*  be  the  centre  of  the 
entrance  pupil.  Placing  an  adequate  optical  system  at  Q*  on  the  optical  axis  with 
Qi,0  as  its  focal  length  /,  we  can  fix  the  following  relations  between  Eo>f?o  and 
E',G',  the  principal  variables  in  front  and  behind  the  optical  system, 

E'g=-J 

G'«=0,  E'k=A'IAo  (1-28) 

where  /=!//.  These  expressions  satisfy  the  relation  G'oE^ K—G'EE'e= A^j At  of 
(1-19),  and  mean  that  G'o  and  S' a  are  constants  and  the  rays  from  a  point  at 
infinity  (£^o=0,  Go=l)  concentrate  to  a  point  Q  through  this  optical  system  as  an 
aberration-free  image,  and  that  a  point  (£o=l.  Go=0)  on  the  entrance  pupil  forms 
its  image  on  the  same  plane  because  of  G'*=0. 

Also  it  is  seen  that  a  point  Q*(£o=li  JGt=A'IAt—\)  on  the  spherical  surface 
that  passes  through  Q*  and  has  its  centre  at  Q  remains  unchanged  because  of  G'=Go 
and  £'£o  as  derived  from  (1-9). 

We  shall  call  such  an  optical  system  the  ideal  optical  system.  A  general  optical 
system  can  be  transformed  into  a  standard  optical  system  by  placing  the  ideal 
optical  system  in  front  of  it,  whereby  our  problem  can  be  treated  in  this  unified 
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way,  and  from  this  standpoint  various  calculations  will  be  dealt  with  hereafter. 

Besides,  relations  that  are  required  in  later  sections  are  derived  by  using 
(1-22)  as 

D'=E'*mDo-2JE'eBo+P{Ho-C*).  (1-29) 

Here  H'=Ho,  and  it  is  obvious  that  the  relations  1~4  stated  in  Section  1-6 

hold  as  they  are.  Hence,  an  incident  ray  in  a  general  optical  system  can  clearly 
be  specified  by  subsidiary  variables. 

Transforming  the  last  formula  of  (1-29)  by  using  S' b= A' jAo  of  (1-28),  we 
obtain 

E'^e-2BJE'  b+P{Ho-C^)=\.  (1-30) 

Solution  of  this  equation  for  S' e  is  given  by 

E'e =Vl-pHoA^o+JBo .  (1-31) 

Of  cause,  for  a  paraxial  ray,  E'e  becomes  1. 

10.  Linearity  of  Aberration  Coefficients 

/ 

It  is  shown  that  aberration  coefficients  of  ordinary  optical  system  can  be  obtained  * 

from  the  aberration  coefficients,  which  is  of  the  same  order,  of  the  standard  optical 
system  with  linear  relations. 

Let  us  use  principal  and  fundamental  variables  in  this  discussion. 

In  the  standard  optical  system,  the  principal  variables  of  the  image  point  are 
expressed  as 

X'G=X'Go+(X'Gaa+X'atb+X'acC) 

-\-iiX'oatO*-\-JCGabOb-\-’ '  ‘  ‘  (1-32) 

where  X'go  is  the  value  of  X'g  for  a  paraxial  ray,  therefore  the  differential  coeffi¬ 
cients  of  X'a  are  regarded  as  the  aberration  coefficients  of  X'a. 

From  the  calculation  about  pupil,  J?'*  is  obtained  which  is  expressed  as 

X'  E=X'  B0'\~{X'  BaO-^Xf  ElJ>-\-X'  BcC) 

-\-{\X’ EafiO}-\-X' EabOb-\-’ .  (1-33) 

In  an  ordinary  optical  system,  of  which  the  principal  variables  are  G  and  E, 

1C =X' gG+X'  eE,  and  the  following  expression  is  obtained  from  the  above  formulas 

X'=X'a+{X'aa+X',;b+X’cc) 
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X’a  =  X’GaG-\-X>BaE,  X* ^=X' X’ e„E,- ■  •  • 
X'a»  =  i(-X^OoiG  +  X^ Ba,*E),  X' at  =  X'aa»G-\-X'  BabE,‘ 
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and  ' . 

ix'-x'^)*X'a’i'  =  -AX'at+{x'+3x'^)X'a6-2x\{x' 
+x'^)(X',>,+X'ac)+x'\(Zx' 

+x'^)X'te-4x'x'>^X'e»- •  •  • .  (1-38) 

From  (1-34)  and  (1-38),  it  is  seen  that  X'a'.X'v,'",  have  linear  relations  with 
X  Oat  "  *  ■  '  >  X' Eat  ••••and  also  that  X'  a*t,  X  a'6S  '  *  •  •have  linear  relations  with  X'aatt 
••••,A''«..,  ••••. 

Hence,  when  the  aberration  coefficients  with  respect  to  the  image  and  pupil 
in  the  standard  optical  system  are  computed,  the  aberration  coefficients  of  the  image 
in  an  ordinary  optical  system  can  be  obtained. 

These  relations  hold  further  among  quantities  linearly  combined  with  principal 
variables  E'  and  G'. 

Fundamental  variables  are  employed  in  this  discussion,  nevertheless,  subsidiary 
variables  can  be  taken  as  such  quantities  by  the  simple  reason  that  fundamental 
variables  are  more  suitable  to  define  incident  rays  in  an  ordinary  optical  system. 

Result 

Characteristic  variables,  namely  the  subsidiary  and  principal  variables,  are 
introduced,  so  that  the  ray  tracing  and  differetial  calculation  will  be  made  simple 
and  easy. 

Especially,  it  is  noted  that  the  linearity  among  principal  variables  will  simplify 
the  differential  calculation  hereafter. 

Besides,  by  employing  more  variables  than  what  are  necessary  to  define  the 
geometrical  condition  of  a  ray  and  the  image  point,  several  checking  formulas  and 
fundamental  relations  can  be  derived. 
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Abstract 

The  derivation  of  ray  tracing  formulas  which  involve  the  case  of  skew 
ray  and  are  suited  to  differential  calculation  is  described. 

The  standard  optical  system  and  the  expression  of  principal  variables  used 
in  Part  I  are  also  employed  in  this  part. 

It  is  shown  that  the  variation  of  the  direction  of  a  ray  in  the  cases  of 
refraction  and  reflection  at  one  surface  and  transfer  from  one  surface  to  the 
next  surface  can  be  expressed  as  the  rotation  of  co-ordinate  system,  one  axis 
of  which  is  the  ray.  According  to  this  fundamental  operation,  the  ray  tracing 
formulas  as  regards  principal  and  subsidiary  variables  are  derived. 

Besides,  formulas  for  an  aspherical  surface  are  obtained  by  extending  the 
above  method  and  formulas  for  a  plane  surface  are  also  derived. 

It  is  shown  by  the  ray  tracing  formulas  that  the  linear  relations  among 
principal  variables  which  are  verified  in  Part  I  hold  throughout  the  optical 
system. 

Simplified  formulas  for  tracing  typical  rays  such  as  paraxial,  axial  and 
extra  axial  rays,  are  derived. 

Further,  choosing  the  wave  front  aberration  as  a  quantity  suited  to  describe 
the  condition  of  the  image,  we  obtain  a  method  to  compute  this  aberration 
with  the  use  of  the  results  of  ray  tracing. 

1.  Introduction 

In  the  cases  of  refraction  and  transfer,  we  begin  with  the  relation  between  the 
variation  of  a  ray  and  that  of  the  co-ordinate  system  -q,  C  on  the  ray. 


\  (  ^21.  f22,  §23) 


>§'t1.§23)  \ 


f. 


Fig.  II-l 
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The  ray  before  and  after  refraction  is  always  in  the  incident  plane,  namely, 
the  ^-Ti  plane,  and  C-axis  remains  unchanged. 

Let  the  angle  of  rotation  about  C-axis  be  Ai,  the  sign  of  the  rotation  from  f- 
axis  to  i?-axis  being  taken  as  positive,  then  the  relation  between  f ,  7  and  tj  becomes 
as  shown  in  Fig.  II-l  and  the  following  relations  hold  among  the  x-,y-,z-  com¬ 
ponents  of  these  co-ordinates. 

p'ii=(/nCosAi+p2,sinAi,  g'2,=flr2iCosAi— pusinAj,  p'*i=(7ji 

a',2=(;i2Cos  At-t-flr22sin  Aj,  (;'22=ff22Cos  Ai— {jf^sin  Ai,  p'j2=(/s2 
g'»=Qa  cos  M+g23  sin  M,  g'2s=g2a  cos  M—gu  sin  Af,  (7'jj=pss  (H-l) 

Next,  for  transfer,  the  direction  of  a  ray  does  not  change  and  the  vzuiation  of 
co-ordinates  are  expressed  with  the  rotation  of  ly-and  C-axis  about  ?-axis. 


’Z'  (  ^21,  fz2,  ^23) 
♦  ^  ?22,  ^23) 


If  (p  is  the  angle  of  rotation,  the  sign  of  rotation  from  y-axis  towards  C-axis 
being  taken  as  positive,  the  relations  between  y,  C  and  r/',  C'  are  shown  in  Fig.  II-2. 

The  relations  among  x-,y-,z-  (Components  of  these  co-ordinates  are  expressed 
as  follows 

gn=g'iu  ff2i=p'2iCos<p-g',isin9;,  gsi=g'si  cos^c-(;'2i  sin^o 

gii=g'i2,  g22=g'22Cosv’— g's2siny),  gs2=g'MCosy)— g'22sin<(> 
gij=g'is.  g2s=g'2»  cos  ^-g'a,  sin  gs3=g'js  cos  (p-g'2a  sin  <P .  (II-2) 

Ray  tracing  formulas  will  be  derived  with  the  use  of  the  fundamental  relations 
derived  in  this  section. 


2.  Formulas  for  Spherical  Surfase 

First,  several  equations  concerning  the  refraction  of  a  ray  at  a  spherical  surface 
are  given  below. 
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Fig.  II-3 

In  Fig.  II-3,  let  the  intersection  between  the  ray  and  the  spherical  surface  be 
R,  and  put  RO=r’  The  sign  of  r  is  the  same  as  that  of  x-component  of  RO. 
Let  the  indices  of  refraction  of  spaces  in  front  and  behind  the  spherical  surface  be 
n  and  n'  respectively,  and  put 

p=nr,  p'=n'r  (11-3) 

where  r  is  a  quantity  which  has  no  concern  with  refractive  index.  Let  OP=h, 
OF=h'  and  the  angles  between  the  normal  at  R  and  the  rays  before  and  after 
refraction  be  i  and  i'  respectively.  Then  we  have 

^sini=A,  p' sin  i'=h',  h'=h. 

Ai=i'— I  is  the  change  of  the  direction  of  the  ray,  and  coincides  with  that 
in  (II-l). 

The  following  quantities  . 

L=coshi,  M=sinhi/h  ’  , 

are  of  trigonometrical  functions,  but  can  be  expressed  more  conveniently  as  follows. 
If  we  put  RP=q  and  RP'=<f  in  Fig.  II-3,  we  have 

q*=p*-H,  q'*=p'*-H 

pp'L=qq'+H.  pp'M=q-q' .  (IM) 

In  the  case  of  reflection,  we  take  «'  =  — n.  The  following  relation  hold  betwen 
L  and  M 

L*+M^H=1.  (U-5) 

From  these  formulas,  it  is  seen  that  L  and  M  are  the  functions  of  H  only. 
Using  the  above  results,  we  derive  the  relations  between  the  principal  or  subsidiary 
variables  concerning  the  ray  before  and  after  refraction,  '  /v/ 
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Combining  (II-l)  with  h  emd  h',  and  using  the  quantities  L  and  M  mentioned 
above,  we  obtain 

g'i\=Lgn+Mhgzi,  h'g'2i=Lhg2i—MHgiu  A'(/'si=Affsi 
g'i2—Lg\z-{-Mhgz2,  h'g,z2=Lhg22 — MHg\z,  h'g'tz^hgzz 
g'i%=Lg\z+Mhgzz,  h'g'zi=Lhgzi—MHgiz,  h'g'zz=hgn.  (II-6) 

From  these  quantities  as  well  as  (1-3)  and  (1-16),  we  have  the  following  trans¬ 
formation  formulas  for  the  principal  and  subsidiary  variables,  If  we  denote  one  of 
the  variables  Eo  and  Go  by  X 

A'=LA-MB.  E'x=LEx+MFx 
B'=LB+MHA,  F'x=LFx-MHEx 

G'x=Gx.  (II-7) 


3.  Formulas  for  Transfer 

In  this  section,  several  formulas  which  hold  for  the  ray  in  passing  from  one 
spherical  surface  to  the  next  are  introduced. 


In  Fig.  II-4,  put  0'0=e,  which  quantity  will  be  made  use  of  hereafter.  Let  the 
distance  between  the  two  spherical  surfaces  measured  on  the  optical  axis  be  d. 
Then  we  have  (the  dash  signifies  the  quantities  for  the  v-th  surface  to  discriminate 
from  those  of  the  v+lth  surface). 

e=n'{d+r—r')  (II-8) 

where  d  has  no  concern  with  refractive  index. 
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Fig.  II-5 

Fig.  II-5  shows  the  projections  of  OP  and  OF,  on  a  plane  perpendicular  to 
OCy,  on  which  the  projection  O  and  O'  are  brought  one  upon  the  other. 

Let  the  difference  in  directions  of  OP  and  OF,  that  is  the  angle  between  ri- 
and  ly'-axes,  be  <p,  which  was  introduced  in  II-§  1.  ly-and  C-components  of  e  are  egzi 
and  epsi  respectively.  Thus  from  Fig.  II-5  it  follows  that 

siny)=efir'si,  -h=egailh',  cosq>={h'—eg'ii)lh={h+egii)lh'.  (II-9) 

This  relations  enable  to  derive  transfer  formulas  for  the  principal  and  sub¬ 
sidiary  variables. 

Combining  (II-2)  with  h  and  h’,  and  using  (11-9)  we  have 
g\i=g'\u  hgii=h'g'ix-e{\-Q'\\\  hgai=h'g'ai 

Pi2=(7'i2,  hg22=h'g't2+eg'iig'ii,  hga2=h'g',2—eg'ia 

g\a=g'n,  hg2a=h'g'23+etfMifix,  hgaa=h'g'ai+eg'x2.  (II-IO) 

Then  using  (1-3)  and  (1-16),  we  obtain  the  following  transfer  formulas  for  the 
principal  and  subsidiary  variables. 

A=A',  Ex=E'x 

D=\-A\  Fx^Fx+eE'xA' 

B=B'+eD',  Gx=G'x+eFx 

H=H'+2eB'+^iy.  (II-ll) 

If  calculations  for  X=Eo  and  X=Go  are  carried  out,  subsidiary  variables  can 
be  computed  from  (1-22)  with  better  accuracy  than  from  the  above  ray  tracing 
calculation. 
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4.  An  Aspherical  Surface 

To  obtain  the  ray  tracing  formulas  for  an  aspherical  surface,  we  begin  with 
the  investigation  of  foundamental  properties  of  an  aspherical  surface  on  the  basis 
of  properties  of  a  spherical  surface. 


Fig.  11-6 


In  Fig.  II-6  which  shows  an  aspherical  surface,  let  the  intersection  between 
the  optical  axis  and  the  normal  to  the  aspherical  surface  at  point  /?  be  O,  and  let 
the  length  of  the  perpendicular  RS  be  y  which  is  drawn  to  the  optical  axis  from 
R.  Put  0'0=p,  PO—p  and  0'S=x. 

Then 

dx 

P^-V'={p-x)*,  {p-xy-^=y.  (11-13) 

Denoting  quantities  for  a  spherical  surface  with  suffices  O  in  this  section,  we 
have  in  this  case  p=p=po  and 

Po'-y'MPo-XoY .  (11-14) 

The  expansion  of  Xo  in  termsi  of  Co={ylpo)*  is  given  by 

XolPo=  - •  (11-15) 

The  term  of  Co  determines  the  property  of  the  paraxial  ray,  but  the  terms  of 
Co*.  Co*. - have  no  connection  with  them. 

Regarding  these  terms  as  characterizing  the  aspherical  surface,  we  express  the 
relation  between  x  and  c  in  this  case  as 

xlpo=  -^c+-^(l-hr,)c*+-^(l+r*)c»+^(H-r,)c*+  •  •  •  •  (11-16) 

where  ri,  r2,***‘may  be  regarded  as  representing  the  degree  of  deviation  from  the 
sphere.  From  this  formula 
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I  \  dx 
K  y  dy 

becomes 

—  =  l+-^(l+t’i)c+-^(H-f2)c*+-^(H-rjV:*+-  •  •  • 

and  (11-13)  is  expressed  as  q—x=k,  p*—y^=ic^. 

If  we  put 

p/Po=l  +  PiC+AiC*+PsC*+ - 

PlPo=l+PiC+p2<:'+p^*+- •  •  •  (11-17) 

then  the  above  relations  give  the  coefficienents  of  the  formulas  (11-17)  as  follows 
Pi  =  —iTu  8p2=5r,+2r,*— 3r2 

16ps=  —6  r,*+2ri»+7r2+6rir2— 5rj 

Pi  =  -iTi,  8/>2=6r,+2ri*-3r2 
16/>s=  — 3ti— 6ti*— 2ri*+9r2+6rir2— Sr* .  (11-18) 

p  and  p  are  computed  easily  and  the  ray  tracing  is  performed  by  regarding  p 
and  p  as  optical  constants  if  c  is  known. 

Now,  c  and  Co  are  expressed  with  the  subsidiary  variables  by  taking  the  origin 
at  the  centre  of  the  sphere  of  radius  po  as 

>l*/>o*c=(x-p)*Z)+2(x-p)B+//-C* 

^Wco=(x-/>o)*Z)+2(x-/»o)B+//-C* 

and  by  subtraction,  we  obtain 

C=Co+iT,c*B'+ •  •  •  •  (11-19) 

where  B'  represents  the  subsidiary  variable  B  when  the  origin  is  taken  at  O'. 
The  second  term  is  an  extremely  small  quantity  in  general. 

If  SO=ST—OT  in  Fig.  II-6  is  expressed  by  subsidiary  variables, 

poW=Po*-(qA+Bf  (11-20) 

is  obtained.  Therefore,  c  can  be  computed  from  (11-20)  and  (11-19). 

5.  Formulas  in  the  Case  when  the  Surface  is  a  Plane 

When  a  plane  surface  is  considered,  the  center,  on  which  the  origin  of  co¬ 
ordinates  is  taken  in  the  case  of  a  spherical  surface,  becomes  indeterminate. 
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In  Fig.  II-7,  let  the  intersection  between  the'plane  and  x-axis  be  the  origin  O, 
and  that  between  the  ray  and  the  plane  be  R. 

P  and  F  are  dealt  with  similarly  to  the  case  of  the  spherical  surface,  in  other 
words  I,  r,  g,  q',  h  and  h'  are  determined  in  the  same  way  as  before.  Let  the  angle 
between  OR  and  the  line  of  intersection  of  the  incident  plane  with  the  plane  be  0, 
and  the  angles  that  OR  makes  with  the  incident  and  refracted  rays  be  u  and  u' 
respectively.  Then  we  obtain 

i4=cosi,  A'=cosi',  D=sin*i,  Z>'=sin*i' 

-qA=B,  -q'A'=B' 

cos  M = COS  ^  sin » ,  cos  «' = cos  O'  sin  i' .  (11-21) 

On  the  other  hand,  considering  OR  as  a  quantity  with  no  regard  to  refractive 
index,  we  have 

n.OR‘Cosu=q,  n'.OR-cosu'=q' 
which  lead  to  ^  ' 

q'lq=n'  cos  u'/n  cos  u=n'  sin  i' In  sin  i=l 
n*(A'*+<7'*)=«'W+9*).  (11-20) 

From  the  relations  obtained  above,  formulas  for  tracing  principal  and  subsidiary 
variables  will  be  derived. 

Put 


n'ln=K,  e=A'IA . 


(11-23) 


Then  from  ,(11-21)  and  (11-22)  we  obtain  for  the  relations  of  D  to  O'  and  B 
to  B’ 


k*D'  =  D,  B'=eB. 


(11-24) 
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Furthermore,  (11-22)  gives  H’ + B'*! A'*=k\H+ A*)  which  becomes  by  the 
relation  of  (1-4)  as 

(11-25) 

We  shall  consider  next  the  formulas  for  the  principal  variables  on  Fig.  II-8, 


In  this  case  QQ'  is  parallel  to  x-axis.  The  expression  (I-ll)  concerning  E 
still  holds  and  moreover  EIE'=n'ln.  Therefore  the  following  relations  hold  for  Q 
and  O' 

nHs'+q')=uHs+q) 
n'x'/cos  r  =  n'*x*/cos  i . 

It  is  seen  from  the  above  that  the  relations  for  the  principal  variables  are 
kE'x  =  Ex, 

rF'x=>ciFx-a-<c*)BExlA 

G'x=k€Gx.  (11-26) 

6.  Paraxial  Ray  and  Axial  Ray 

By  the  use  of  the  ray  tracing  formulas  derived  above,  formulas  for  a  paraxial 
ray  and  an  axial  ray  are  obtained. 

In  the  case  of  a  paraxial  ray,  subsidiary  and  fundamental  variables  are  zero, 
which  make  g=p,  and  therefore  L—\  in  (II-4). 

From  (1-26)  and  (1-25)  it  follows  that  D=0,  Gx—Fx-  The  ray  tracing  formulas 
(II-7)  then  become 

M=llp'-lip,  E'x=Ex+MGx.  G'x=Gx.  (11-27) 

As  for  transfer,  from  (II-ll)  the  ray  tracing  formulas  are  given  by 
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E^x  =  E^ X,  G+i=G'x-\-€E' X  •  (11-28) 

Next  the  method  of  calculation  for  an  axial  ray  is  considered. 


In  this  case,  Bo=Do=0  and,  H  and  D  can  be  easily  obtained  from  (1-22)  as 
H=G(^Ho,  D=Eo*Ho.  (11-29) 

Then  L,  Af  and  A  are  calculated  from  (II-4).  As  is  shown  in  Fig.  II-9,  the 
object  point  Q  and  the  image  point  O'  are  on  the  optical  axis,  and  as  explained  in 
Section  1-3,  Fe=GaA.  Thus  from  (II-7),  the  ray  tracing  formulas  are 

E^a=EEa'\~MGaA,  G.^a—Go-\-cE.^a  (11-30) 

where  the  suffix  +  signifies  the  quantities  for  the  next  surface. 

The  difference  between  —GojEa  of  an  axial  ray  and  that  of  a  paraxial  ray 
represents  the  spherical  aberration  and  the  difference  between  —IjEa  of  an  axial 
ray  and  that  of  a  paraxial  ray  represents  the  deviation  from  sine  condition. 

In  the  general  case  involving  skew  rays,  ray  tracing  can  be  performed  with 
only  Ex  and  Gx  excepting  Fx  as' treated  above,  since  Fx  can  be  computed  from 
(1-24),  (1-26)  and  (1-27)  in  case  of  need. 

7.  Extra  Axial  Ray 

We  treat  here  the  method  of  obtaining  the  distortion,  sagittal  and  meridional 
image  surfaces  from  the  ray  tracing  about  an  extra  axial  principal  ray. 

Using  Ee  and  Ge  instead  of  Ea  and  Ga  dealt  with  in  the  previous  section,  we 
perform  calculations  with  respect  to  the  pupil.  That  is,  we  have  instead  of  (11-29) 
and  (11-30) 

H=GhDo,  D=E>eDo 
E+E=LEE-\-MAGEt  G+b=Ge~\'sE^e  • 


(11-31) 
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Fig.  II-IO 


Since  9  in  Fig.  II-IO  can  be  obtained  from  the  relation  A=cos9,  the  point 
O'  where  the  ray  crosses  the  ideal  image  surface  is  determined.  In  other  words, 
the  distortion  is  calculated.  From  the  results,  Ea  and  Ga  are  obtained  by  using 
the  formulas 

AE^t-a— A*Eg-\- MGa,  G+o=Go+c£+o  (11-32) 

which  are  derived  from  (II-7)  and  (II-ll).  Since  Ea  and  Go,  which  give  the  posi¬ 
tion  of  Q,  are  computed  from  the  above  formulas,  and  since  Ee  and  Ge  for  Q*  are 
already  known.  O'  and  therefore  distortion  are  obtained  with  better  accuracy. 

Calculations  about  Ea  and  Go  correspond  to  those  of  the  sagittal  image  surface 
as  mentioned  in  Section  1-3. 

The  calculation  of  the  meridional  image  can  be  performed  with  formulas  derived 
below  by  using  the  results  obtained  above. 


0  Jt  0-^ 

Fig-ill-ll 


In  Fig.  11-11,  let  Q  and  O'  be  the  meridional  object  and  image  points  of  this 
surface  respectively  and,  PQ  and  PQ'  be  the  lengths  with  no  regard  to  refractive 
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index,  then  we  have 

1/m'  •  /V  =  1/m  •  1  ’q-Hq' 

=  1^-00*  cos  6' . 

If  these  quantities  are  expressed  as  n' •P'Q'  =—G'JE\,  n-PQ=—G^lE^,  Vq'— 
llq=M^,  then  the  above  formulas  become 

E' i^=E^-\- M^Gi^,  G'*=G* 

Et^=E't,  Gt^=G\+eA'E'^.  (11-33) 

In  this  case,  there  are  also  linear  relations  among  the  variables  corresponding 
to  the  principal  variables  as  seen  from  the  above  formulas. 

Of  cause,  for  a  paraxial  ray  (11-33)  coincides  with  (11-27)  and  (11-28). 


8.  Wave  Front  Aberration 

We  shall  derive  the  wave  front  aberration  for  the  ideal  image  point  of  rays 
from  an  object  point. 


In  Fig.  11-12,  let  O  be  the  centre  of  the  exist  pupil,  and  Q  be  the  ideal  image 
point  where  the  wave  front  aberration  is  considered. 

Let  O,  be  the  intersection  between  the  reference  sphere  passing  through  O  and 
having  the  radius  OQ=p  and  the  ray  passing  through  the  image  point  O'!.  Let 
the  perpendicular  from  Q  to  the  ray  be  QQ"i.  If  QQ"i  is  a  small  quantity,  we 
may  assume  that  Q"i  and  Q  are  at  equal  distances  from  Oi. 
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Next,  take  a  neighbouring  ray  Q'2Q"2  shown  in  Fig.  11-13  and  denote 
Q^'i  =  5.  'i^2Q"2  =  S  +  dS. 

On  the  ray  Q'2^'2,  let  the  point  on  the  same  wave  front  as  that  at  O'!  be  at 
a  distance  ds'  from  Q't.  Similarly,  let  the  point  on  the  same  wave  front  as  that 
at  Q"i  be  at  a  distance  dW  from  Q"t.  Then  we  have 

dW=ds-ds' .  (11-34) 

According  to  this  basic  idea,  we  shall  find  the  fundamental  formulas  expressing 
the  differential  dW  concerning  the  wave  front  on  Cartesian  co-ordinates.  We  shall 
obtain  s  from  Q{x,y,z),  Q’{x',v',2f)  and  the  direction  cosines  of  the  ray  gi,gi,g3 
in  Fig.  13.  Since  s  in  the  component  of  Q'lQ  in  the  direction  of  Q'iQ"t  we  have 

s={x-x')gi+(y-if)g2+{z-3^)gi .  (11-35) 

If  Q  is  assumed  unchanged,  differentiation  of  this  expression  gives 

ds=ix-x')dgi-\-(.y-V')dg2-\-{z—z!)dgt 
-igidx' ->rgidy'  +gtdz:) . 

Next,  d^  is  expressed  as  the  variation  of  co-ordinates  of  and  Q't  as 
ds'  =  -igidx'  +g2dy'  +gsd^) . 

Consequently,  when  these  two  equations  are  combined,  dW  given  by  (11-34) 
becomes 

dW=ix-x')dgi+{y-y')dgi+{2-^)dg, .  (11-36) 

Now,  we  proceed  to  express  the  fundamental  formula  obtained  above  in  terms 
of  variables  so  far  employed. 
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Representing  Eg  and  Go  as  principal  variables  for  and  £oo  and  Gao  for  Q, 
we  have 

X=— Gao!  Ego,  I/  =  W^oo,  z=2oIEgo 

3^ =— Gal  Eg,  y'=yolEG,  z'=2oIEg.  (11-37) 

Let  the  differences  be 

X— x'=A(Gg/£g) 
y-l/'  =  -A(l/£o)l/o 

z-z!=-^(MEg)zo.  .  (11-38) 

Then,  expressing  gi,g2  and  g^  in  terms  of  the  principal  variables,  we  modify 
(1-14)  with  (11-37),  namely 

gi=A,  togi=EEyo—E(ni*o,  toOi=EEZo—EGZto.  (11-39) 

Using  these  quantities  in  the  above  fundamental  formula  (11-36)  and  assuming 
that  yo  and  Zo  are  constants,  we  have 

dW=^{G/E)dA+Hl|E)dp  (11-40) 

where  P=(£g'6o— £«Oo)//o.  When  the  expressions  of  (1-21)  are  used,  it  is  seen  that 
P  coincides  with  Pa  in  (1-26). 

For  an  axial  ray,  p=Q,  and  (11-40)  becomes 

dW=KGIE)dA.  (11-41) 

Practical  value  of  W  can  be  computed  from  the  following  integral  by  using 

the  results  of  ray  tracing. 

I 

W=^^{<GIE)dA+^l^(XIE)dP  (11-42) 

or  for  an  axial  ray 

W=^^{GIE)dA .  (IM3) 

Result 

Ray  tracing  formuljis  as  regards  principal  and  subsidiary  variables  and  the 
formula  of  wave  front  aberration  are  derived  with  which  calculation  for  skew  rays 
can  be  dealt  with  simply. 
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It  is  noted  that  calculation  of  imaginary  rays  which  have  no  physical  meanings, 
for  example  in  the  case  of  Do=Bo=0,  Ho<l  and  Do=Ho=0  Bo^O,  is  possible  in 
the  same  way  as  mentioned  above. 

Therefoire,  complicated  character  of  imagery  can  be  analized  to  a  certain  extent 
by  means  of  the  calculation  of  imaginary  rays. 
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Abstract 

The  differential  calculation  applicable  for  obtaining  aberration  coefficients, 
which  have  been  derived  by  other  methods,  is  described,  • 

As  mentioned  in  Part  I.  subsidiary  variables  are  adequate  quantities  to 
define  the  incident  ray,  to  separate  aberration  terms  and  to  differentiate  the 
ray  tracing  formulas. 

Therefore,  differentiating  the  ray  tracing  formulas  by  the  subsidiary  vari¬ 
ables  in  the  object  space,  we  can  derive  the  differential  formulas  of  ray  tracing 
and  use  them  to  calculate  the  derivatives  of  principal  variables  successively 
from  surface  to  surface. 

These  derivatives  contain  terms  of  higher  order  of  the  subsidiary  variables 
for  the  incident  ray,  but  by  using  the  paraxial  quantities  in  the  above  deriva¬ 
tives,  we  can  obtain  the  paraxial  differentiation  that  is  the  differentiation 
in  the  region  of  Seidel. 

However,  in  this  treatise,  formulas  different  from  the  above  ray  tracing 
formulas  are  used. 

By  successive  differentiations,  the  2nd  order  coefficients  and  the  3rd  order 
coefficients,  the  latter  of  which  have  scarcely  been  dealt  with  in  practice,  are 
obtained. 

When  such  differential  calculation  is  extended  to  the  wave  front  aberration, 
the  differential  coefficient  thus  obtained  is  the  required  aberration  coefficient. 

Further,  several  methods  to  separate  and  compute  the  aberration  terms  by 
performing  the  calculation  on  imaginary  rays,  are  briefly  referred  to. 

\i  Introduction”*’ 

General  consideration  about  differential  calculation  is  given  in  this  section. 

A  quantity,  for  exzunple  the  principal  variable  in  the  image  space  represented 
by  X',  can  be  expanded  in  terms  of  the  subsidiary  variables  in  the  object  space  Ha, 
Bn  and  Do  £is 

X' = X'«Do) 

-{-{iX'n*  hbHoBo-\~  •  •  •  •)+ '  *  •  •  (lll-l) 

where  X'o  represents  the  paraxial  value  of  X'. 

1)  F.  Wachendorf :  Optik,  5  (1949)  80. 

2)  H.  A.  Buchdahl ;  Optical  Aberration  Coefficients,  (Oxford  University  Press.  London,  1954). 
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If  all  the  differential  coefficients  are  equal  to  zero,  then  X’  =  X'o  and  the  image 
becomes  ideal.  Hence  it  is  reasonable  to  regard  these  coefficients  as  aberration 
coefficients  of  X'. 

The  differentiation  of  X'  by  one  of  the  subsidiary  variables,  for  example 
Ho,  gives 

dX'/dHo=X'HMX'HfHo=X'HBBo+X'i,DDo)+--  •  • .  (III-2) 

The  derivative  obtained  from  differentiation  of  ray  tracing  formulas  is  of  the 
form  given  above,  and  let  the  differentiation  of  this  form  be  called  the  differentia¬ 
tion  of  ray  tracing  formulas. 

The  derivative  dX'IdHo  involves  various  terms  of  higher  order  of  subsidiary 
variables,  but  if  paraxial  values  are  introduced  (that  means  Ho=Bo=Do=0),  then 

dX'ldHo=X'H  (III-3) 

is  obtained.  Such  a  differentiation  is  called  the  paraxial  differentiation.  This 
treatment  corresponds  to  that  in  the  Seidel’s  region. 

2.  Diflferentiation  of  Ray  Tracing  Formulas 

Derivatives  of  principal  and  subsidiary  variables  by  one  of  the  subsidiary 
variables  in  the  object  space  represented  by  a  are  derived  in  the  case  of  a  spherical 
surface. 

Since  L  and  M  in  (II-4)  are  functions  of  only  H,  differentiation  of  these  by 
H  gives 

2qqH=-l,  2q'q'H=—l 
PffMH=^qH-q' H,  PP'LH=(fqH+qq'H  +  \ 

Nh=M^MhH  (III-4) 

where  N=MH.  Then  from  (II-7)  and  (U-ll)  the  derivatives  of  subsidiary  variables 
by  a  are  obtained  as 

A\=LA.-MB.HLhA-MhB)H. 

H.=LB.+NA.->r{LHB+NHA)H. 

and 

A^.=A'.,  O'. - 2A'A\ 

B^.=H.+eD', 

H^.=H'.+2eB'.+^D'..  (III-5) 

By  the  use  of  these  formulas,  it  is  possible  to  compute  the  derivatives  A*„ 
B^m.D^m  and  H*.  with  respect  to  the  next  surface- from  the  derivatives  A,,B,  and 
Hm  with  respect  to  this  surface. 
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Next  from  (11-32)  the  derivatives  of  principal  variables  are 

AE-t-Xm'^  A.^Exm-¥MG Xm-\-\E^^Xm 

G^Xm  =  Gxm-\-€E^Xm 

where 

[E^\xm= A*mEx-\-A,E*.x-\-MHHmGx .  (III-6) 

When  the  derivatives  Exm  and  Gxm  for  the  incident  ray  to  this  surface  are 
zero,  only  remains  which  may  be  therefore  called  the  peculiar  term  to  this 

surface. 

3.  Paraxial  Differentiation 

As  mentioned  in  the  abstract,  when  paraxial  values  are  taken  in  the  ordinary 
differentiation,  the  differentiation  is  limited  in  the  region  of  the  1st  order.  How¬ 
ever,  for  convenience  sake,  we  employ  formulas  different  from  those  used  in  the 
previous  section. 

First  we  take  up  formulas  for  H  and  D  from  (1-22) 

H  =G*  f.Dq-¥2G  F.GaBo-\~G^(AHn — C*)-t-C* 
Z)=£*j?Z)o+2£K£ofio+£^fi(^o — C*)  • 

Since  it  is  seen  from  (II-4)  and  (1-3)  that  M  is  a  function  of  only  H  and  A 
is  a  function  of  only  D,  they  may  be  expanded  as 

o 

pp'M,=p-p',  2pp'M„=llp'-\lp 

From  (1-22)  the  derivatives  of  subsidiary  variables  are  given  by 
G^eDoa  -f  2G£Go-Bo«  ^■G*oZ)o« 

Da=E^EDoa-\-2EEEaBoa-\-E^aDo<, .  (I1I-8) 

In  practice,  there  remains  only  one  term  among  the  terms  involving  Do,,  Bo, 
and  Do,,  for  example,  if  a=Ho,  then  Bo,=Do,=0. 

Next,  by  differentiating  (III-7),  Af,  and  A,  are  obtained,  namely 

M,=MhH„  A,=  -\D,.  (III-9) 

Differentiating  (11-32)  with  ..the  use  of  these  relations  and  bringing  the  terms 
peculiar  to  this  surface  together  as  done  in  the  previous  section,  we  have 
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\E*\xm=A-^mEx~AmE^X-\-MmGx  •  (III-IO) 

This  is  what  corresponds  to  the  so-called  Seidel’s  term  peculiar  to  this  surface. 
Then  from  (11-32),  the  derivatives  of  principal  variables  are  derived  as  follows 

E^Xm=Exm-\-  MGxm-\-{E.^\xm 

G*Xm  =  Gxm-\-eE.t.Xa  •  (IIHl) 

The  derivatives  of  principal  variables  in  image  space  may  be  computed  by 
applying  these  linear  type  formulas  from  surface  to  surface. 

According  to  this  reason,  if  the  paraxial  derivative  of  A’  by  E'  is  represented 
by  X'e’,  X' xm  can  be  expressed  as  the  following  formula,  which  is  the  sum  of 
derivatives  relating  to  all  surfaces. 

X'x.=  i:X'x’[EAx..  (111-12) 

It  is  seen  from  the  above  relations  that  X'  xm  can  be  computed  only  with 
paraxial  values. 


4.  Second  and  Third  Order  Paraxial  Differentiations 

We  proceed  to  investigate  the  paraxial  differentiation  of  the  2nd  order  by 
applying  the  results  of  and  the  method  employed  for  the  1st  order  differentiation. 

Differentiation  of  (III-8)  by  which  represents  one  of  the  subsidiary  variables 
Ha,  Bo  and  Do,  gives 

\H mf—G  BmDof¥GamBo^ 

+Go2(G««fio#+ GomHof) + i  ( 1 — G*a)C*mf 
iDmf=EBj^(EBmDof-i-EamBof) 

+Ea^(EB^Bof-\-EaaHof) — hE}oC*mf.  (III-13)  * 

The  notation  £  means,  for  example,  Z;Gjr.£)o^=GK.A^+Gjs^Do..  Then  from 
(III-9) 

A.f-=-kD.f-\D.Df  (III-14) 

are  derived.  From  (III-ll)  the  derivatives  of  principal  variables  become 

E.^Xmf=Exmf-^MGxmf-\-[E^.]xmf 

G.^Xm^=Gxm^-\■eE^Xm^ 

where 

{Ef\xaf={A.t^mfEx  —  A.fE^X  —  Aa^E^X-\-M,fGx) 

+  2](A  +  ,£jr^  —  AaE^Xf-^-MmGxf)  •  (III-15) 


Expressed  in  sums,  it  becomes 
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X’x.f=ZX'AEAx.^.  (III-16) 

Further,  we  derive  derivatives  of  the  3rd  order  employing  the  results  of  the 
1st  and  2nd  order  differentiations. 

Differentiation  of  (III-13)  by  r  gives 

=  G  EnfDo,  +  GomfiBor) 

+G«2](Gi,^D+Go«f//or) 

+  2  (G  e  ,G  ^  Dor  +  (xErnGn^Bor  +  GomGofHtr) 

—  2]  GoGo«C*^r 
and 

J  D,pr  =  EE^{EEt$Dor  +  EcmpBor) 

+  Ea^iEEmpBor  +  EompHor) 

+  ^iEEmEEpDoT-\'EE»EapBoj-\-  Ed.EupHor) 

-ZEaEo.Op,.  (III-17) 


Then  from  (III-14) 

M.pr-^MnH.pr  +  MH^ZHMpr  +  MmHMpHr 

A.p,=  -hD.p,-\T.D.Dpr-lD.DpDr  (III-18) 

are  obtained.  By  differentiating  (III-15),  formulas  of  ray  tracing  type  is  derived  as 

E^,XmPJ=ExmPr^\■MGxmPJ-\■\E^.\xm^t 

G^XmPT  =  GxmPr'\-6E^XmPr  (III-19) 

where  the  term  peculiar  to  this  surface  is 

[^+]  XmPr  =  (A^mPjEx  —  AmPjE^.X-\-Map,Gx) 

+  l^(,A^,pExj  —  AmpE^Xj-\-MmlGxr) 

+  YLiA^..Expr~AmE^xpr'^-MmGxpr)  •  (III-20) 

t 

Expression  in  sum  form  becomes 

X'x.pr=ZX'E'[EAx.pr.  (III-21) 


5.  Differentiation  in  the  Case  when  the 
Surface  is  a  Plane 

Formulas  for  computation  in  the  case  of  a  plane  is  summerized  in  this  section. 
Using  (11-24)  and  (11-23)  as  the  basic  formulas,  we  introduce  first  the  deriva¬ 
tives  of  subsidiary  variables,  namely 

x^D'  =  D,  A^^=l-iy,  e=A'IA. 

Differentiation  of  these  by  a  gives  O'.,  A',  and  «,  as  follows 
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2A’A'.=  -iy.,  e.le=A'JA\-A.IA.  (III-22) 

Differentiation  of  the  above  gives 

kE^.x=Ex,  G*jr=«Gjr+e£+jr  (III-23) 

which  are  obtained  by  combining  (11-26)  and  (II-ll),  and  denoting  the  term 
peculiar  to  this  surface  as 

[G+]jr«=<'*,G.r  (III-24) 

we  obtain 

^E.^Xm  =  Exm 

G^xm=rtGxm-\-eE^xm-^[G4.]xm  •  (III-25) 


For  paraxial  differentiations,  if  we  differentiate  (11-23)  by  a,  $  and  r  employing 
the  derivatives  of  A  given  in  Section  3  and  Section  4,  the  derivatives  of  «  are 
obtained  as 

9m“A-^m~Am9  ea0~A^mf  Amf 

tmfr~  A.t.mfr  —  Amh —  ^**Afj —  ^fmfAr  •  (III-26) 

Results  of  successive  differentiations  with  the  use  of  the  above  relations  are 
given  below. 

By  differentiation  of  the  1st  order 

[G^]xm  =  etmGx,  eE^Xm  —  Exm 

G^Xm  =  KGxm-¥eE^Xm-¥\Gf\xm  •  (III-27) 

By  differentiation  of  the  2nd  order 

\Gf\xmf=rt,ffjxAf7^*mGxf  eE*Xm^=Exmt  , 

G.t^Xmf  =  KGxmf-\-eE^Xmf+[G*.]xmf  .  (III-28) 

By  differentiation  of  the  3rd  order 

[G.t]xmfr=Kt,fjGx-\'e^t,fGxt-\-K^tmGxfr 

rE*Xmff=ExmfT 

G*Xmfr  —  '^rmfT'^^^*^‘tr'^lG-¥]xmfT  •  (III-29) 

6.  Differentiation  for  Ideal  Optical  System 

The  differential  expression  for  the  ideal  optical  system  is  sought  to  establish 
the  calculation  method  to  be  dealt  with  as  of  the  case  of  the  standard  optical 
system. 

The  derivatives  of  principal  variables  by  a  are  given  from  (1-28)  as 
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G'o.=£'o.=G'*.=0  (III-30) 

and  from  (1-30)  as 

2JiE'  e-JB,)E'  b.=2JE'  eBo.-JKH^-0.)  .  (III-31) 

The  derivatives  of  subsidiary  variables  by  a  are  obtained  from  (1-22)  as  men¬ 
tioned  before  by  using  the  above  values. 

For  paraxial  differentiation,  (11-30)  also  holds  and  the  successive  derivatives  of 
E'  E  are 

E'  Em  =  A,  —  E'  EA^m  , 

E' Emf=Aaf  —  E' eAq^P — ^E' E*Aof  , 

E'  EmPr~  AmPr~E'  EAomff —  J^E'  EmAoPr —  YiB'  EmpAfij  .  (III-32) 


7.  CoeflScients  of  Wave  Fr<Hit  Aberration 

The  wave  front  aberration  coefficients  of  the  1st,  2nd  and  3rd  orders  will  be 
obtained  by  the  use  of  the  fundamental  formulas  derived  in  Section  II-8  and  the 
results  of  paraxial  differentiation. 

Differentials  of  (11-40)  can  be  expressed  with  respect  to  Ho  and  Bo  can  be  ex¬ 
pressed  as  follows : 

WV/= A(Go/£o)/4//-l-A(l/£o)P// , 
WB=t^{GalEo)AB+KVEo)PB.  (III-33) 

Values  of  Wh  and  Wb  can  be  computed  from  the  differentiation  of  ray  tracing, 
hence  W  can  be  obtained  from  the  following  integral  with  respect  to  Ho  and  Bo : 

W=  J  WHdHo+ 1  WBdBo .  (III-34) 

Since  Q  in  Fig.  11-12  is  the  ideal  image  point,  if  we  use  paraxial  values  in 
(III-33)  we  have  Wh=  Wb=0.  Hence  there  is  no  1st  order  differential  coefficients 
of  W,  and  by  differentiation  of  (III-33)  with  respect  to  a,  we  have 

WHm=(GalEa)aAH 

WBm={GalEa\AB-VO-IEG)mPB  (III-35) 

in  which  six  coefficients  relating  to  the  three  cases  of  a=Ho,  =  Bo  and  =Do  appear 
but  not  the  coefficient  Wot.  Whb  is  obtained  from  either  one  of  the  above  formulas 
and  the  relation  shown  below  holds,  therefore  there  are  five  independent  coefficients 
corresponding  to  the  five  terms  of  Seidel’s  aberrations. 

From  (III-9),  (III-8)  and  (1-26)  A,  P  and  therefore  the  wave  front  aberrations 
are  expressed  with  the  principal  variables.  In  this  way  the  relations 
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2WHt=—EaGaH,  W^hb=— 

^Whb— — E(^aB,  IVb>= — EgG'oB 
2Whd= — EdGoDy  Wbd— — EeG'od 

W„=0  (m-36) 

are  derived  where  G'  represents  the  principal  variable  G  when  the  centre  of  the 
pupil  is  taken  as  the  origin. 

Next,  differentiation  of  (III-35)  with  respect  to  $  gives  the  following  aberration 
coefficients  of  the  2nd  order : 

(Go!  Ea)mfA  H 

■\-Yt{GalEa)mAHf+YtO-l^d).PHf  , 

WBmf  =  {Gal  Ea)mfA  B  +  (XIEa)mfPB 

+  ^{Ga/ Ea)mA Bf-\-  J10-IEa),PBf .  (III-37) 

Similarly,  the  3rd  order  aberration  coefficients  are  obtained  by  differentiation 
of  (III-37)  with  respect  to  r  thus 

Wnmfr^  (Gn/ Ea)mfTA  u 

+  IHGalEa)mtA  ^OIEG)mfPHr 

+  YtiGalEa)tAH^1-¥  HO-IEo)mPHh  y 

WBmfr  =  {Go!  Ea)mfrA  B + 

+  ^{Gal  Ea)mtA  Br"^  ^0-1  Ea)mfPBT 
+  Jl(.Ga/Ea)y,ABfT-{-  Yf\IPo)mPBh  •  (III-38) 

')  • 

8.  Classification  and  Separation  of  Aberration  Terms 

By  cleissifying  the  wave  front  aberration  terms  and  aberration  terms  of 
principal  variables  that  will  be  dealt  with  below  and  by  using  the  calculation 
method  given  heretofore,  let  us  consider  the  method  of  separating  these  terms. 

Here  (III-l)  is  again  considered,  namely 

+{.^X' h*Ho^-¥\X' h*bH^Bo-¥"  ’ 

If  all  the  differential  coefficients  are  equal  to  zero,  then  X'^X\  and  -the 
image  becomes  ideal  as  mentioned  before,  and  we  can  regard  these  coefficients 
as  aberration  coefficients  which  can  be  classified  as  arranged  in  the  following 
diagram.  •  ,;  •>  '  •; 
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In  the  following,  we  consider  the  application  of  the  previously  introduced 
calculation  to  these  aberration  terms. 

First  if  X+  is  the  value  of  X  when  Bo=Do=0  and  //o>0,  we  have 

X'+=X'o-\-X' HHo-\-iX' H*Ho^'^iX' ‘  ‘  .  (III-39) 

Similarly,  as  to  the  imaginary  ray  for  which  Ho  is  negative  and  has  the  same 
absolute  value  as  used  above, 

X'-=X'o-X'HHo+iX'HtHo^-iX'„,Ho^+ . .  •  •  (III-40) 

is  derived  in  which  X'  is  represented  by  X'-.  Addition  and  subtraction  of  these 
two  formulas  give 

^(,X\-X'.)=X^HHo^-^X’H^Ho^+^X'HtHo'‘+--  •  (III-41) 

respectively,  from  which  we  see  that  odd  and  even  coefficients  can  be  separated, 
and  further,  if  we  compute  X'  for  several  values  of  H,  coefficients  of  higher  order 
can  be  sepcUBted. 

Extending  this  method  of  calculation  to  an  imaginary  ray  such  as  Z>o=0, 
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Bo^O,  Ho^O,  we  can  separate  the  coefficients  involving  Ho  and  Bq. 

If  the  derivative  of  X'  with  respect  to  Bo  represented  by  dXjdBo  is  computed 
by  the  method  of  ordinary  differentiation,  the  following  expression  is  obtained : 

dXldBo=X'B-^X„BHo+  \X'„2bHo^+-  •  •  •  .  (III-42) 

X'b.X'hb.Xh^b, - can  be  separated  by  the  method  given  above. 

By  the  repetition  of  such  treatment,  each  group  of  terms  lying  on  the  dotted 
line  in  Diagram  III-l  can  be  separated. 

It  is  seen  that,  in  the  case  of  a  paraxial  differentiation,  coefficients  X'h,X'b 
and  X  D  on  the  1st  zone  are  determined  by  the  differential  calculation  of  the  1st 

order,  X'ih,X'hb, - on  the  2nd  zone  by  the  differential  calculation  of  the  2nd 

order  and  so  on. 

As  stated  before,  there  exists  neither  the  coefficient  of  the  wave  front  aberration 
concerning  only  Do  nor  the  coefficient  of  the  1st  order  corresponding  to  the  term 
on  the  1st  zone  in  Diagram  III-l.  Accordingly  there  appear  five  aberration  coeffi¬ 
cients  of  the  1st  order,  nine  coefficients  of  the  2nd  order  and  fourteen  coefficients 
of  the  3rd  order  corresponding  to  the  terms  on  the  2nd,  3rd  and  fourth  zones 
respectively. 


Results 

Two  kinds  of  differential  calculation  methods  are  derived. 

The  first  that  is  the  differentiation  of  ray  tracing  formulas  is  available  to 
obtain  detailed  information  about  the  complicated  aberrations  of  an  optical  system. 

The  second  that  is  the  paraxial  differentiation  is  used  to  obtain  aberration 
coefficients  of  the  quantity  suited  to  describe  the  imagery  directly  as  the  differential 
coefficients.  ' 

It  seems  difficult  to  use  such  a  differential  method  with  respect  to  the  ordinary 
ray  tracing  formulas  which  are  generally  being  used. 
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Abstract 

For  calculation  of  correction  on  an  optical  system,  various  methods  of 
calculation,  which  depend  on  the  nature  of  the  correction,  are  being  used. 

In  this  treatise,  a  synthesized  differential  method  of  calculation  is  con¬ 
sidered  which  is  applicable  to  the  cases  of  an  aspherical  surface,  an  eccentric 
optical  system,  displacements  of  object  plane  and  entrance  pupil,  and  so  on. 

The  case  of  differentiation  of  ray  tracing  formula  and  the  case  of  paraxial 
differentiation  are  dealt  with  both  involving  skew  rays. 

When  any  part  of  the  optical  system  deviates  by  a  small  amount,  the 
principal  and  subsidiary  variables  in  this  space  change  and  the  principal  vari¬ 
ables  in  image  space  also  change  by  the  amount  determined  by  the  partial 
differential  relations  among  these  variables. 

Although,  for  this  partial  differential  calculation,  the  method,  by  which  the 
differential  coefficients  of  principal  variables  in  image  space  with  respect  to 
variables  in  any  space  successively  from  image  to  object  space  is  computed,  is  » 

commonly  being  used,  here  we  compute  the  coefficients  directly  using  the  * 

results  of  differential  calculations  obtained  in  Part  III. 

After  this  treatment,  calculation  formula  of  differential  correction  is  derived 
in  the  case  in  which  surface  elements  vary,  or  in  other  words,  optical  con¬ 
stants  vary. 

Next  for  an  eccentric  optical  system,  in  which  the  center  of  any  surface 
deviates  from  the  optical  axis,  and  for  an  aspherical  surface,  the  surface  which 
was  originally  spherical,  formulas  of  corrections  are  derived. 

Further,  the  correction  formulas  for  an  optical  system,  in  which  the  posi¬ 
tions  of  object  plane  and  entrance  pupil  are  displaced,  are  obtained. 

Lastly,  formulas  of  differential  correction  for  the  wave  front  aberration 
concerning  the  above  various  variations  are  expressed  in  general  form. 

1.  Introduction 

Partial  differential  calculations  have  been  employed  for  corrections  of  an  optical 
system  in  general.”*’ 

First,  the  principle  of  the  method  of  obtaining  the  relation  between  partial 
differential  calculation  and  correcting  calculation,  is  considered  by  employing  the 
notations  hitherto  used. 

1)  F.  D.  Cruickshank:  Proc.  Phys.  Soc.  57,  (1945)  350. 

2)  Tatsuro  Suzuki :  Text  Book  of  Kogaku  Konwa  Kai  (1953). 


40 


A  Differential  Method  of  Optical  Calculation.  IV 


41 


Derivatives  of  X'  by  X  is  expressed  as  follows  in  terms  of  the  partial  deriva¬ 
tives  of  X'  with  respect  to  E’  and  G'  which  are  the  principal  variables  of  the 
emerging  ray  of  the  vth  surface. 

X'x  —  X'fE'x+X'a’G'x.  (IV-1) 

When  a  quantity  K  of  a  part  of  the  optical  system  is  varied  by  a  small  amount 
dV,  an  emerging  ray  and  finally  X'x  are  changed  by  the  amount  defined  by  the 
following  relation 

X'  xv  =  X'  pE'  xv-i-X'a’G'  xv+X'  pvE'  x  +  X'o'vG'  X 

where 

X'  PV  =  X*  PH'II'v-\-X'  e<b'B'v\-X'  PD'IX  V 

X'a  v=X'a'H'H'v-\-X'a'B’B'v\-X'o'D’D'v  (IV-2) 

which  are  derived  from  the  differentiation  of  (IV-1)  by  V. 

Therefore,  in  ordinary  differentiation,  X'xv  can  be  computed  when  the  deriva¬ 
tives  of  X'  and  those  with  respect  to  V  are  obtained  in  the  above  formulas. 

For  a  paraxial  ray,  H'v.B'v  and  IXv  in  (IV-2)  become  zero,  then 

X' XV  =  X' B'E' XY X' O'G' XV  (IV-3) 

is  derived. 

As  for  paraxial  differentiation  to  obtaine  X'xv.,  the  derivative  of  X’xm  by  V, 
we  differentiate  first  (IV-l)  by  a,  thus 

X'  Xm  =  X'  pE'  Xm  +  X'o'G'  Xm-\-X'  p,E'  X-\-X'a‘<,G'  X 

where 

X'  Pm  =  X'  PH’D'  ,-\-X'  E'B'B'  m-\-X'  PD’D' t 

X’a’.=X'o'H>H'.+X'a  pB'.+X’a'D>D'. .  (lV-4) 

Then  the  differentiation  of  the  above  by  V  gives  X'x.v,  namely 

X'  XmV  =  X'  pE'  XmV+X'a'G'  XmvhX'  P  mE' XV Xo’ aG' XV 
+  X'  PmvE' xAX'a'mvG' X 

where 

X'  PmV  =  X'  E’H'H'  mV-\-X'  E'B'B'  mvAX'  E'D'IX  ,V 

X' O'mV  —  X' mv¥X' iq-bB' mvVX' a>D>iy mV  •  (IV-5) 

These  formulas  coincide  with  what  are  derived  from  the  differentiation  of 
(IV-2)  by  a. 

The  calculation  method  of  the  above  partial  differential  coefficients  of  X'  will 
be  dealt  with  in  the  following  section. 


2.  Partial  Differential  Calculation 

The  method,  by  which  the  partial  differential  calculations  of  principal  variables 
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are  carried  out  by  the  use  of  the  results  of  differential  calculations  of  ray  tracing, 
will  be  considered  first. 

Using  the  relation  of  (1-19),  namely 

G(ioEi>o—GboE{io=AIAo  (IV-6) 


and  solving  the  equation  of  (IV-1)  for  X=Eo  and  X=Go  as 
X  BQ^X'nGEO'i'X'EEEO 

X'm  =  X'nGoo + X'  eE/io 

we  have 

X' «  =  (G«oX'  eo-GeoXoo) 
Xa=^  (EEoX'ao-En^X'Eo). 


(IV-7) 


(IV-8) 


The  right  terms  of  these  formulas  can  be  computed  from  the  previous  dif¬ 
ferential  calculation.  For  a  paraxial  ray,  A  =  Ao=l  and  the  above  formulas  become 


X'  e  =  GooX'  EQ  —  GEoXno 

Xa=EEoX’oo-EooX'Eo.  (IV-9) 

The  partial  derivatives  by  subsidiary  variables  are  obtained  similarly.  To  do 
this,  (IV-1)  is  differentiated  and  the  derivative  of  X'  by  X  and  a  is  expressed 
through  E,  G,  H,  B  and  D  as 

X'  X. =X’eEx.+X\/;x.+X’  e.Ex + X'c..Gx 

where 

X'  Et  =  X’  EHHm  +  X'  EliB,-{-X  EnDm 

Xo,=X'ohH.+X’obB.+X'odD,  .  (IV-10) 


X' E.  and  X'am,  which  will  be  used  later  on,  are  obtained  as  follows. 

$ 

If  we  put 

{X']x.=X  x.-X'  eEx.-X'oGx.  (IV-11) 

this  is  a  quantity  obtainable  by  the  calculations  mentioned  before,  with  which 
(IV-10)  for  the  cases  of  X=Eo  and  X=Go  are  expressed  as  follows 
X'  EmEEO'hX'naGEO=  [X']eo« 

X'omGao+X' E,iEao=[X]ao<i  • 

By  combining  these  two  formulas  and  using  the  relation  of  (1-19),  we  have 
X’e.  =  ^  (Gao[X']Eo.-GEo[X']„o.) 

X'«.=  ~{Ee,XX'\o^-Eoo{X'\eo.)  .  (lV-12) 
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Next,  for  a=Ho,Bo  and  Do,  the  formula  of  X' g.,  which  is  the  first  formula  of 
(IV-10),  becomes 

X'  EhHho'^’X'  EbBu^'^X'  edDho  =  X'  KHt 

X  EBHBO-k-X'  EBBBn-\-X'  EdDbo—X'  EBO 

X  EhHdo-\-X' bbBdo-\-X'  eI>Ddo  —  X'  EDO  •  (IV-13) 

X'eh.Xeb  and  Xed  (similarly  XaH,XoB  and  Xod)  are  obtained  by  solving 
these  simultaneous  equations.  For  paraxial  differentiation,  the  following  relation 
derived  from  (1-22)  is  used  to  solve  the  above  equations: 

Hho  Bho  Dho 

Hbo  Bbc  D«„  =1.  (IV-14) 

Hdo  Hdo  Hdo 

The  derivatives  by  E',G',H',B'  and  O'  concerning  the  emerging  ray  are 
obtained  in  the  same  was  as  above  by  using  E'  =  E+,  G'=G,  0^  =  0^  and  H'  —  H 
and  calculating  B',  from  (1-22)  similarly  to  the  case  of  //,  and  D.. 

3.  Preliminary  Calculation 

If  the  distance  from  the  origin  of  one  surface  to  the  origin  of  the  next  surface 
is  varied,  the  relation  of  the  incident  ray  to  the  latter  surface  will  be  changed. 

The  derivative  of  X x=X' o-^G^x'\'X^ e*E.^x  of  (IV-1)  by  e  is 

X'  x»  =  X'u*G^.Xt  +  X  O'^G^Xt  +  X'xH*  lI*»-\-X'  XB^-B*e 

X'  X  D^D*t 

where 

X'  XH*  =  X'o*H*G^X-‘tX'  E*H*E^X  ' 

X'xB^  =  X'  a*B*G^X  +  X'  E*B*E+x 

X' XD*  =  Xa*D*'\-G^X-¥X  e*d*E*x  •  (IV-15) 

Since  the  derivatives  of  X'  in  the  right-hand  side  of  (IV-15)  are  obtained  by 
the  previous  method,  X xt  can  be  computed  if  the  derivatives  by  e  are  known. 

Differentiation  of  (II-ll)  by  e  gives 

G^X*  =  E.t.X,  E^.Xe  =  0 

//♦,=2B+,  B^,  =  D^,  D^e=0.  (IV-16) 

Using  these  quantities,  we  obtain  the  required  relation 
,  X' xt^X'u*E^x-\-2X' xh*B^-{-X xb*D^  .  (lV-17) 

For  a  paraxial  ray,  Bi.=D^=0  and  the  above  becomes 
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X'x,=X'o*E^x .  (IV-18) 

In  the  case  of  paraxial  differentiation,  we  obtain  X'  xm  and  differentiate  it  by 
e,  but  the  result  is  the  same  when  (IV-17)  is  differentiated  by  a  as  explained  in 
Introduction,  therefore  we  have 

X’ x,.=X'a.E^x.+X'o^.E^B+2X' xh^B^.+X’ xb.D^.  .  (IV-19) 

Next,  the  variations  of  the  optical  constants  concerning  the  surface  are  ex¬ 
pressed  as  the  differentials  of  L  and  M,  but  there  exists  the  relation 
of  (II-5),  the  differentiation  of  which  gives 

LLm=-MH.  '  (IV-20) 

If  the  formulas  of  (II-7)  are  regarded  as  expressing  the  relations  between  the 
variables  about  the  rays  before  and  after  passing  the  surface  and  if  they  are  dif¬ 
ferentiated  by  L,  M  and  the  above  quantity  Lm  is  used,  we  have 

G'xm=Q  AE'xm  =  ExA'm+G'x 

//'*=0,  B'm=H'A’IL.  A'm  =  -B'IL,  IXh^-IA'A’m.  (IV-21) 

The  derivative  of  X'x  by  M  is  expressed  through  the  variables  about  the 
emerging  ray  from  this  surface  as 

X'  XM==X'g'G'  XM-^X'  E'E'  XM-\-X'  XH'H'  M 
+X'xb'B'm+X'  XD>  IXm 

where 

X'xH'  =  X'E-H'E'x  +  X'a'H'G'x 

X'xB'=X'e'B'E'x  +  X'o'B>G’x 

X' xD'  =  X' e'D’E' x-\-X'a’D>G' X  (IV-22) 

and  from  (IV-21)  and  (IV-22),  the  required  differential  formula  is  obtained  as 
X  XM — X'  e’{ExA'  m~\'G'  x)IA+{X'  xb'H'  -^■2X'  xd’B')A'  IL .  (IV-23) 

For  a  paraxial  ray,  since  the  subsidiary  variables  do  not  take  part  in  the 
differentiation, 

X'xm=X'eG'x.  (IV-24) 

The  paraxial  differentiation  is  the  differentiation  of  (IV-23)  by  o,  namely 

X'  Xt,M=X'  E'{G'  Xt  —  B'  aEx—G'  xAm)-\-X'  B'mG'  X-{-X'  XB'B’  . 

^2X'  xd'D'  . . 

In  this  case  the  variation  of  the  surface  involves  the  variation  of  Mh  in 
(III-9),  and  therefore  the  differentiation  by  Mu  gives  the  following  derivative 
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X' XmMH  —  X' fE' XmtlH  =  X' E-G' xH' m  .  (lV-26) 

Relations  derived  in  this  section  will  be  widely  applied  hereafter  to  the  cal¬ 
culation  of  corrections. 


4.  Differential  Corrections  for  Surface  Elements 

The  radius  of  a  surface,  the  refractive  index  of  a  space  and  differential  varia¬ 
tions  of  e,  M  and  Mh,  which  are  the  quantities  derived  from  the  surface  elements, 
will  be  considered  for  obtaining  differential  variation  of  X' x  with  res:ect  to  the 
above  quantities  by  using  the  results  of  the  previous  section.  X'  x  by  the  above 
surface  elements  will  be  obtained  from  the  results  in  the  previous  section. 

First,  differentiations  of  (II-8)  by  a,n',r  and  r*  give 

ed=^n',  en'=eln',  e,=  —n',  ej^=n' .  (IV-27) 

Differentials  of  (II-4)  by  r,  n  and  namely  the  derivatives  of  Af  by  these 
quantities,  are  given  by 

PP'Mj = npjq — n’p'l<t — 2nn'rM 

pp'Mn=riPtq-Mp').  pp'Mn’ - riPiQ-Mp).  (IV-28) 

For  a  paraxial  ray,  we  put  p=q,  P'=q'  in  the  above  formulas,  which  become 
Mr=-Mlr,  Mn=Vnp,  Mn’  =  -lln'p'.  (IV-29) 

Differentiation  of  (III-7)  gives  the  derivatives  of  Mh  expressed  as 

2pp'MHr=—jM,  2npt/MHn=-2lp-ltp' 

2n'pp'MHn'=llp-2lp' .  (IV-30) 


p-/  V'  iz-t-  I 
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In  Fig.  IV-1,  the  relations  between  the  quantities  in  front  and  behined  the 
vth  surface,  for  which  variations  are  considered,  are  shown.  The  variation  of  d 
has  effect  upon  e,  that  of  r  upon  e-,  e+,  M  and  Mh,  that  of  n'  upon  e,  M.  M+, 
Mh  and  M*h^.  Therefore  in  an  ordinary  case,  using  the  quantities  of  (IV-17), 
(lV-23),  (IV-27)  and  (IV-28),  we  obtain 

X'  Yd  =  X'  XrC,,  X'  Xf  =  {X'  XeCT*)-~X'  XtCf-h  X'  X  ¥ 

X'xn’  =  X'  XeSn’-^X'  XMMn’  +  iX'  X¥Mn)+  (IV-31) 

where  (  )_  and  (  )+  are  the  quantities  concerning  the  v-lth  and  the  v+lth 

surface  respectively. 

For  paraxial  differentiation,  similarly  to  the  above,  we  obtain  the  differential 
correction  formulas  for  X'x  as 

X'  Xot  —  X'  XarCd  , 

X’ XiMj  —  iX' Xarej  +  )--\-X' XafCj  +  X' XmmMt  +  X' XaMHMuj  , 

X'  X  on'  =  X’  x.^n'  +  X'  Xo¥Mn'MX'  Xo¥Mn)* 

X' Xo¥nMHn'-\-{X' XoMH^Hn’)*  .  (lV-32) 


5.  Elccentric  Optical  System 

From  the  results  of  calculations  concerning  the  partial  derivatives  in  Section 
IV-2,  formulas  for  calculation  of  derivatives  of  X'x  are  derived  for  the  case  in 
which  the  center  of  one  of  the  spherical  surfaces  in  the  optical  system  is  displaced 
from  the  optical  axis. 

Consider  a  co-ordinate  system  with  the  centre  at  the  centre  O  of  vth  surface 
of  the  optical  system  described  in  Section  I-l.  For  a  general  treatment  of  an 
eccentric  optical  system,  it  is  sufficient  if  we  consider  that  the  fundamental  plane 
varies  and  that  the  displacement  of  the  centre  of  the  vth  surface  is  only  in  y- 
direction. 

The  position  of  the  object  point  Q(x,y,z)  relative  to  the  vth  surface  as  well 
as  the  position  of  the  point  Qn(Xn,y*,Zn)  relative  to  the  entrance  pupil  remain 
unaltered. 

However  if  Ay  is  the  deviation  of  the  centre,  y  and  y*  become  yolEa+^y  and 
ynolEe+^y  respectively.  It  is  seen  from  these  expressions  that  Eg  and  Eb  may 
be  considered  to  be  constant.  In  other  words,  the  principal  variables  do  not  change 
with  the  deviation  of  the  centre  from  the  optical  axis. 
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Among  the  subsidiary  variables,  D  depends  only  on  the  direction  of  the  ray 
and  is  unchanged,  but  //  and  B  vary  because  the  perpendicular  OP  from  O  in 
Fig.  IV-2  changes. 

The  variations  of  H  and  B  are  calculated  by  the  use  of  (1-14)  as 

#W=(x2+a)  (x**+c)-(xx*+6)* , 

—  x*a+(xx*+6)— xc , 

a=y^+z^,  b=yyit+zzi^,  c=v**+x** ,  (IV-33) 

and  the  altered  values  ‘of  variables  in  image  space  are  computed  with  ray  tracing 
formula. 

For  differential  correction,  expressing  the  deviation  as  a  differential  dy=dy^ 
and  using  x,Xi^,t  which  are  unchanged,  we  obtain  by  differentiating  the  above  ex¬ 
pressions  with  respect  to  y  or  y* 

= (x^* -  (XX* -f  6)  (y y*)-f  (X* -f  a)y  * 
/*B,=(x-x*)(y-y*).  (IV-34) 

In  the  case  of  ordinary  differentiation  these  quantities  are  computed  from  the 
relations  in  Section  1-5.  Then  for  paraxial  differentiation,  we  put  a=b=c=0  and 
X— x*=f,  and  we  have  for  //*  and  By 

iff/,=xy,-x*y,  tBy=y-y^.  (IV-35) 

They  may  also  be  expressed  in  terms  of  principal  variables  as 

iHy=—GEyo+Gnyto,  By=Eay*o-EByo.  (IV-36) 

X'xyt  the  differential  variation  of  X’x  caused  by  the  deriation,  is  expressed 

through  H  and  B  as 


X' Xy=X' xnHy  +  X' xaBy  . 


(IV-37) 
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X'xH  and  X'xB  are  the  quantities  obtained  from  the  partial  differential  cal¬ 
culations  stated  before. 

For  paraxial  differentiation,  using  (IV-36)  in  the  above  formula,  we  obtain 
X'  xt=(—2GEX'  xh-^EbX'  xB)yo'^i^OoX'  xH—EaX'  xB)ymo  •  (IV-38) 

Calculation  of  these  coefficients  with  respect  to  Vo  and  |/*o  gives  the  quantities 
corresponding  to  the  aberration  coefficients  caused  by  the  deviation  of  the  centre. 

The  above  is  the  case  when  the  optical  system  behind  the  vth  surface  is  dis¬ 
placed  by  the  amount  dy^=dy.  If  the  variation  caused  by  the  displacement  rfi/v+i 
=  —dy  of  the  optical  system  behind  the  v4-lth  surface  is  added  to  the  above,  the 
variation  oi  X'  x  caused  by  the  displacement  of  only  one  spherical  surface  is  given. 

Thus 

^'x,=X'x,.-X'xy.*x .  (IV-39) 

Besides,  the  variation  of  y'  and  y\  in  image  space  are  expressed  by  the  com¬ 
bination  of  quantities  corresponding  to  the  magnification  and  focal  length  of  the 
paraxial  ray 

y'y=(E'ya-E\^  ia)IE'a 

y'^,*={E\E-E\^xE)IE'E.  (IV-40) 

6.  Aspherical  Surface  Modification 

With  the  use  of  the  results  of  partial  differentiation  and  fundamental  relations 
for  an  aspherical  surface  derived  in  Section  II-4,  differential  formulas  in  the  case 
in  which  any  sphere  constituting  the  optical  system  is  subjected  to  an  infinitesimal 
change  and  becomes  aspherical,  but  retains  the  rotational  symmetry  about  the 
optical  axis,  will  be  considered. 
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The  discussion  will  be  made  on  the  practical  stand-point. 

First  from  (11-17)  and  (11-18),  neglecting  terms  of  c*,c®, — ,  we  have 

P=P=Po — (IV-41) 

2po 

Next  we  proceed  to  express  y*  in  terms  of  variables  that  have  been  employed. 
By  taking  c=Co  of  (11-19)  and  using  (11-20), 

y^=Po'-(qA+Br  (IV-12) 

is  obtained.  When  transformed  by  cf=p^—H  and  when  only  the  terms  of  the  1st 
order  is  taken,  the  above  becomes 

y*=i/-2po5+po*7) . 

By  comparing  this  with  (II-ll),  we  see  that  y*  corresponds  to  H  in  the  case 
in  which  the  origin  is  displaced  by  a  distance  — Po  and  brought  to  O'.  That  is  to 
say,  taking  O'F^  in  the  above  figure  as  H\,  we  have 

y2=//''> .  (IV-43) 

As  is  shown  in  (IV-41),  p  and  p  are  subjected  to  the  same  change  in  this 
aspherical  surface.  This  corresponds  to  the  variation  of  the  radius  in  the  case  of 
sphere.  Therefore,  — H’^dr  is  the  quantity  that  corresponds  to  the  variation 

^ftPQ 

dr  of  the  radius  in  the  case  of  sphere. 

Using  this  relation,  we  obtain  the  variation  of  X'  against  r  as  . 

X'xr  =  -  -^X'xrH'o.  (IV-44) 

^tiPo 

For  paraxial  differentiation,  it  becomes 

X'x.r - -^X'xrH'o.  .  (IV-45) 

This  is  the  differential  expression  corresponding  to  the  variation  of  Seidel’s 
aberration  term. 


7.  Displacements  of  Object  Plane  and  Elntrance  Pupil 

The  differential  relations  in  the  case,  in  which  the  axial  point  of  the  object 
plane  and  the  entrance  pupil  of  the  ideal  optical  system  treated  in  Section  1-9  are 
displaced,  will  be  considered. 

The  relations  among  variables  in  the  ideal  optical  system  depend  only  on  /  as 
can  be  seen  from  (1-28)  and  (1-29). 
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First,  by  differentiating  these,  we  obtain 

G'aj—0,  E'aj=l,  G' Ej=0 

{E' b-JBo)E' ej  =  E' eBo-F{Ho-0)  .  (IV-46) 

For  subsidiary  variables,  differentiation  of  (1-28)  with  the  use  of  (III-9)  and 
(1-28)  gives  the  following  expressions 

B'j=E'ejBo-Ho+C^ 

\D'j=-E'eE'ejA^^.  (IV-47) 

In  the  case  of  paraxial  differentiation,  the  differentiations  of  (1-28)  and  (1-29) 
by  a  lead  to 

G'aj=G'Ej=E'EJ=0,  E'aj  =  l, 

G'omJ  =  E'otiJ  =G' EaJ  =  0,  E' EtJ  =  Boa  —JH^n  , 

B\j=-Ho.,  iiy.j  =  -E'E.j.  (IV-48) 

For  a=Do,  the  derivatives  by  a  become  zero. 

The  displacements  of  the  object  plane  and  the  entrance  pupil  are  now  con¬ 
sidered  on  Fig.  IV-4  which  corresponds  to  Fig.  1-6  of  the  ease  of  the  ideal  optical 
system. 


Since  the  displacement  of  the  object  plane  corresponds  to  the  variation  of 
/(==!//).  differentiation  of  X’ x=X' e'E' x+X' a'G' x  by  /  gives 

X'  xj  =  X'  a’G'  xj-\-X'  e’E'  xj-\-X'  xh'H'j 

+X'xB'B'j-\-X'xD’iyj .  (IV-49) 

Using  (IV-46)  and  (IV-47)  in  this  formula,  we  obtain  X'xj. 

For  paraxial  differentiation,  differentiation  of  (IV-49)  by  a  gives 
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X'j,J  =  X'  a’G'xmJ  +  X'g'E'xmJ 

-\-X'xH>H'mJ-]-X'xB’B'.J-\-X'xD'D',J 
-\rX'  XH’jB'  a-\-X'  XB’jB'  m  +  X’  XD’jly  m 

where 

X'xH'J  =  X'  E'wE'xmJ  +  X'  O'H'G'xJ 
X'xB’J  =  X'e’B’B'xJ  +  X'qib'G'xJ 

X’xD'J  =  X'e'D'  E'xj-¥Xa'D'G'xj .  (IV-50) 

8.  Differential  Correction  of  Wave  Front  Aberration 

If  the  quantity  concerning  differential  correction,  for  example,  the  surface 
element,  eccentric  deviation,  coefficient  of  an  aspherical  surface,  is  expressed  by  V, 
the  derivative  of  X'x  by  V  that  is  X'xr  and  the  paraxial  derivative  X'xmv  are 
obtained  by  the  methods  mentioned  previously. 

For  ordinary  differentiation,  the  following  formula  of  (11-42) 

W=^MGolEa)dA  +  ^^a/Eo)dP 

is  used.  For  differentiation  of  ^{GalEa)  and  A(l/£o),  differentiation  of  a  paraxial 
ray  is  used  at  the  same  time. 

With  the  differential  notations  hitherto  used,  the  differentiation  of  W  is  ex¬ 
pressed  as 

H'V=  ^^{Gal  Ea)vdA-\-  Ea)vdP 

-\-^MGalEa)dAv+^MllEa)dPr  (IV-51) 

from  which  Wr  can  be  computed. 

For  paraxial  differentiation,  we  have  from  the  formula  of  (III-35) 

WHm  =  {GalEa)mAH  , 

lVB.=(Go/Eo).AB+a/Eo).PB  (IV-52) 

and  their  derivatives  by  V  are  given  by 

WHmV  =  {GalEa)mvAH-\-{GalEa)mAHY  , 

Wfl.  K = {Gal  Ea)m  vA  B-\-{Ga/ Eg)*  A  bv 

+ (1/ Ea)mrP B + (XIEa)mPBV  .  (IV-53) 

For  eccentric  deviation  in  the  case  of  paraxial  differentiation,  X! xy  has  been 
obtained  already,  hence  we  begin  with  the  following  formula  of  (III-33) 
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WH=^{GolEa)AH  > 

W^B=A(Go/£o)i4B+A(l/£o)Pi»  •  (lV-54) 

The  differentiations  of  these  by  y  gives 
WHy  =  (.Gal Ea)yA  H  , 

WBy=(GalEa)tAB+{XIEfi)tPB  •  (IV-55) 


Result 

Differential  correcting  formulas  in  the  case  of  ray  tracing  for  all  kinds  of  rays 
involving  skew  rays  and  imaginary  rays,  are  derived. 

Besides,  correcting  formulas  of  Seidel’s  aberrations  which  can  be  computed  by  ' 
using  only  the  values  obtained  from  the  calculation  for  a  paraxial  ray,  are  derived. 

It  is  shown  that  the  partial  differential  coefficients  required  for  this  calculation 
can  be  computed  directly  from  the  differential  calculation  mentioned  in  Part  III. 

Therefore,  the  differential  calculation  given  in  this  treatise  becomes  most  effec¬ 
tive  when  applied  to  the  calculation  of  aberration  and  that  of  differential  correction. 
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Abstract 

Focal  surface  and  reciprocal  linear  dispersion  of  Ebert  spectrograph  are 
described.  The  focal  surface  varies  with  the  position  of  the  plane  grating. 

When  the  distance  between  the  spherical  mirror  and  the  plane  grating  is 

/?(1— -^y  ),  where  R  is  the  radius  of  the  mirror,  the  focal  surface  is  a  plane. 

When  the  distance  is  ^/?,  the  direction  of  the  dispersion  is  perpendicular  to 
the  grating  grooves.  The  reciprocal  linear  dispersion  varies  with  the  inclina¬ 
tion  of  the  grating  and  also  by  the  position  on  the  focal  surface. 

1.  Introduction 

Optical  properties  of  Ebert  monochromator  when  it  is  used  in  conjunction  with 
a  spectrograph  are  discussed.  An  attempt  was  made  to  clarify  by  calculation  the 
relation  of  the  position  of  grating  to  the  curvature  of  focal  curve  and  the  depen¬ 
dency  of  dispersion  on  diffraction  angle. 

2.  Tangential  and  sagittal  focal  curves 

In  Fig.  1,  O  is  the  center  of  the  spherical  mirror  and  G  is  an  arbitrary  point 


Fig.  1. 
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on  the  line  MO.  The  origin  of  the  rectangular  coordinate  is  at  O.  X-axis  coincides 
with  MO  and  T-axis  is  perpendicular  to  it.  Let  GP,G'F  be  two  parallel  rays  in 
XF-plane  intersecting  at  F  after  reflection  at  the  mirror.  Then  the  coordinates 
of  F{X,  Y)  are 


where  R  is  the  radius  of  the  mirror,  g  is  the  X-coordinate  of  G,  and  x  is  the  angle 
PGM  in  the  figure. 

In  Fig.  2,  G'P'  and  G”F'  are  two  rays,  both  parallel  to  PG,  lying  in  a  plane 
perpendicular  to  XF-plane.  After  reflection  they  intersect  at  F'  which  is  in  XT- 
plane.  The  coordinates  of  F(X',  Y')  are 

^'  =  -f+|-(l-^)sin=z+-.-.  (3) 

Y'=  y  sinz+-j sin»x+-*-' .  (4) 


3.^  [GratingX^^uAtion 

In  Fig.  3,  TZ-plane  is  the  grating  surface  and  the  direction  of  the  grooves  is 
parallel  with  Z-axis.  Let  a  beam  of  parallel  rays,  the  direction  cosines  of  which  are 

— cosC  cos/i,  — sinC  cosm,  sinyu 

be  incident  on  the  grating  and  let  direction  cosines  of  the  diffracted  rays  be 
cos  TQ  cos  V,  sin  n  cos  v,  sin  v. 
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The  optical  path  difference  between  two  rays  diffracted  from  two  adjacent 
grooves  must  be  a  multiple  of  the  wavelength.  Hence,  the  difference  is 

d  sinC  cos fi+d  siniy  cosv=w^  (5) 

where  d  is  the  grating  constant. 

Z 


The  optical  difference  of  two  rays  diffracted  from  the  same  groove  must  be 
zero.  Hence, 


sin/i+sin  v=0 

From  (5)  and  (6).  next  two  relations  are  obtained 

mk 


sin  C+sim7= 


d  cosyi 


n=—v. 


(6) 

(7) 

(8) 


4.  Focal  Curve  of  Spectrograph 

In  Fig.  4,  P  is  a  point  light  source  and  Q  is  its  image.  The  grating  grooves 
are  parallel  with  Z-axis.  Among  the  angles  a,  ip,  ^  and  ^  in  the  figure  there  exist 
two  relations.  When  the  inclination  angle  of  the  grating  normal  is  d. 


sin  (a  -  ») + sin  (^  -  ^)  =  . 

acosy> 

'  9=-^ 

If  a+^>2^,  m  is  positive  and  if  m  is  negative. 


(9) 

(10) 
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Tnagents  of  angles  e  and  e'  which  the  planes  OGP  and  OGQ  make  respectively 
with  Xy-plane  are 


tan  e  = 


tan  (p 


tan  e'  = 


tan^ 

sin^ 


(11) 


Cosines  of  angles  x  and  x'  which  GP  and  GQ'  make  respectively  with  J^-axis  are 
cos  x= cos  a  cosy)  cos  x'=cos  0  cos^ .  (12) 


(9)  and  (10)  show  that  if  the  slit  is  parallel  with  Z-axis,  spectral  lines  appear 
nearly  parallel  with  Z-axis  and  the  spectrum  is  nearly  parallel  with  K-axis. 


V 


In  the  following,  the  position  of  the  slit  is  assumed  to  be  near  X-axis  and  its 
direction  coincide  with  that  of  the  Z-axis.  In  such  a  case  the  rays  which  are 
incident  upon  the  garting  can  be  regarded  as  parallel  rays  and  so  can  be  also  the 
diffracted  rays.  Taking  into  consideration  that  the  direction  of  the  slit  is  parallel 
with  that  of  Z-axis  and  also  that  the  angle  ^  of  the  diffracted  rays  is  small,  we 
obtain  from  (1),  (2),  (11)  and  (12)  the  coordinates  Q{X,  Y,  Z)  of  the  image  point 
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A’=-y+y(l — (sin=^+sin*y5) 

Y=  sin^  cos¥)+|-|-+-^j (sin2^+sin*<p)sin^ 
Z=  - y  sin^p- ^  (sin*  ^4-sin* ip)  sin  if . 


then  ^=“Y' 


(13) 

(14) 

(15) 


This  means  that  the  focal  surface  is  a  plane,  but  as  the  Z-coordinate  is  a  function 
of  the  spectrum  is  curved  in  ZZ-plane. 


that  is,  if  the  grating  is  at  a  distance  -7-  from  the  spherical  mirror,  the  Z-coordi- 

4 

nate  becomes  independent  of  In  this  case  the  spectrum  is  perfectly  parallel  with 
ATT-plane.  This  position  of  grating  is  being  used  in  some  commercial  spectro¬ 
graphs.” 


5.  Dispersion 

3 

When  g  is  equal  to— -7-^ 

4 

.Y=  —1/?—^  /?(sin-  sin-  y) 

>'=-  sin  cos y- 

Z= — ^y?sin  V 

therefore 

rfS=(rfX-  t  sin^j*  f  yCos*y|*/?co8^rf^ 

From  (9)  and  (16)  the  reciprocal  linear  dispersion  is  derived  as 

cos#+sin(!>.-^jl+2-^Jj  . 

The  reciprocal  linear  dispersion  decreases  as  0  increases. 

1)  R.  F.  Jarrel  J.O.S.A.  45  (1955)  259. 


(IS  mR 


(16) 
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As  the  distance  of  the  slit  from  XK-plane  increases,  the  dispersion  varies. 

The  wavelength  4  of  the  rays  that  focus  on  the  line  1^=0  of  the  focal 
surface  is 

(l-ySin*,.).  (17) 

When  ip  is  very  small,  io  may  be  considered  constant  but,  with  the  increase 
of  if,  spectral  lines  begin  to  till  against  XZ-plane. 
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